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STATEMENT  OF  PROBLEM  STUDIED 


Rotational  and  vibrational  collisional  energy  transfer  processes  in  polyatomic  molecules 
have  been  investigated  by  means  of  infrared/millimeter  wave  double  resonance  (IRMMWDR) 
techniques.  Previous  studies  have  generally  been  limited  to  empirical  descriptions  of  a  limited 
subset  of  the  most  predominant  channels  because  of  the  plethora  of  energetically  available 
channels.  An  objective  of  the  present  study  was  to  account  for  all  significant  relaxation  pathways, 
and  to  describe  them  in  terms  of  a  few,  hopefully  physically  meaningful,  parameters.  To  this  end, 
several  isotopic  species  of  CH3F  and  CH3CI  were  studied  in  detail.  More  than  350  IRMMWDR 
time-resolved  transitions  were  measured  in  CH3CI  alone.  These  experimental  observations 
constituted  an  extremely  diverse  set  of  data;  they  included  more  than  40  probe  transitions  with 
various  angular  momentum  quantum  states  J  and  K ,  a  number  of  different  pump  transitions  and 
vibrational  states,  and  a  wide  range  of  pressures.  The  parameters  were  determined  by  fitting  the 
theoretical  model  on  which  they  were  based  to  the  experimental  data  set  by  means  of  a  nonlinear 
least  square?  fit.  All  of  the  rotational  collisional  transitions  could  be  described  using  only  five 
parameters,  and  these  parameters  are  directly  related  to  well  known  physical  parameters. 

Several  features  of  the  experimental  apparatus  were  crucial  to  the  success  of  trie 
experiment.  First,  the  millimeter-wave/submillimeter-wave  probe  is  extremely  sensitive  to  any 
deviations  of  rotational  state  population  distributions  from  ambient,  Boltzmann  distributions.  This 
extreme  sensitivity  allow  ed  the  measurements  of  small  perturbations  in  the  populations  of  states  far 
from  the  pumped  states.  The  sensitivity  also  provided  a  demanding  test  of  the  thermal  nature  of 
several  of  the  theoretically  proposed  rotational  energy  transfer  processes.  Second,  the  essentially 
continuous  tunability  of  the  millimeter-wave/submillimeter-wave  probe  allowed  selection  of  a  wide 
range  of  rotational  transitions.  The  resulting  large  data  set  provided  much  redundancy  in  verifying 
the  theoretical  model  for  the  collisional  energy  transfer.  Additional  details  on  these  and  other 
features  of  the  experimental  set-up  are  described  in  more  depth  in  appendix  A. 

SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

The  most  rapid  collisional  process  is  that  which  changes  the  principal  rotation  quantum 
number  J,  while  not  changing  any  other  quantum  numbers.  Owing  to  the  large  permanent  electric 
dipole  moment,  the  most  rapid  of  these  is  the  dipole  collision  process  whose  selection  rules  are 
AJ  =  0,  ±1.  The  cross-section  for  this  process  can  be  nearly  an  order-of-magnitude  larger  than  the 
gas  kinetic  cross-section.  Although  significantly  slower,  higher  order  rates  AJ  >  1  were  found  to 
make  significant  contributions  to  rotational  energy  transfer  because  of  their  ability  to  affect  large 
changes  in  quantum  numbers  in  a  single  collision.  All  of  the  J-changing  collisions  were  well 
modeled  by  both  the  statistical  power  gap  fitting  law  and  the  infinite  order  sudden  approximation 
with  the  assumption  of  a  power  law  expression  for  the  basis  rates.  Remarkably,  by  simply 
accounting  for  the  difference  in  average  thermal  velocities  between  CH3F  and  CH3CI,  the  same 
values  for  the  pure  AJ  changing  process  could  be  used  for  both  molecules. 


Unlike  with  the  AJ  processes,  significant  differences  were  found  between  the  AK 
processes  of  CH3F  and  CH3CI  (the  projection  of  the  angular  momentum  on  the  molecular  axis  is 
given  by  K).  For  symmetric  top  molecules  with  C3u  symmetry  such  as  CH3F  and  CH3CI, 
collisional  rates  which  change  K  by  other  than  multiples  of  three  are  extraordinarily  slow  since  they 
require  changes  in  the  nuclear  spin,  which  is  very  weakly  coupled  to  molecular  collisions.  The 
earlier  work  on  CH3F  demonstrated  that  rotational  energy  transfer  involving  AK  processes  were 
not  selective  in  l  or  K  aside  from  the  above  requirement  that  AK  =  3n.  In  addition  it  was  found 
that  these  AK  collisions  were  thermal  in  the  sense  that  they  distributed  the  transferred  population 
according  with  a  Boltzmann  distribution  at  a  single,  ambient,  thermal  temperature.  This  was 
reflected  in  the  experimental  observation  that  all  of  the  time-responses  from  the  AK  process  had 
identical  shapes.  In  contrast,  it  was  found  that  the  AK  time-responses  for  CH3CI  varied 
significantly;  exhibiting  differing  degrees  of  a  sharp  feature  which  followed  the  termination  of  the 
pump  pulse.  This  difference  led  to  a  generalization  of  the  model  for  AK  collisions.  We  found  that 
all  of  the  AK  curves  could  be  accurately  modeled  by  simply  allowing  a  time-varying  temperature 
for  die  Boltzmann  distribution.  The  time-varying  temperature  has  the  appealing  features  of  being 
an  exponential  decay  to  ambient  temperature  from  an  initial  temperature  that  is  related  to  the 
difference  between  the  energies  of  the  pumped  state  and  ambient.  This  straightforward  and 
physically  appealing  model  not  only  accounted  for  a  wide  variety  of  AK  data  in  CH3CI,  but  it  also 
satisfactorily  reproduced  the  CH3F  AK  data.  The  experimental  and  theoretical  details  of  these 
results  are  discussed  in  depth  in  appendix  A. 
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Time-resolved  double  resonance  study  of  J-  and  K-changing  rotational 
collisional  processes  in  CH3CI 

Travis  W.  Pape,  Frank  C.  De  Lucia, a)  and  David  D.  Skatrudb) 

Department  of  Physics,  Duke  University,  Box  90305,  Durham,  North  Carolina  27708 
(Received  20  August  1993;  accepted  4  January  1994) 

Time-resolved  double  resonance  spectroscopy  using  infrared  pump  radiation  and  millimeter-wave 
and  submillimeter-wave  probe  radiation  (IRMMDR)  has  been  used  to  study  rotational  energy 
transfer  (RET)  in  CH3C1.  A  collisional  energy  transfer  model  using  only  five  parameters  for  RET 
plus  those  needed  for  vibrational  processes  is  shown  to  accurately  model  350  IRMMDR  time 
responses  for  two  different  pump  states  and  43  probe  transitions  covering  a  wide  range  of  rotational 
states.  Previous  studies  in  this  laboratory  have  revealed  that  J-  and  /(-changing  RET  have  vastly 
different  characters  in  CH3F  [J.  Chem.  Phys.  92,  6480  (1990)].  Both  J-  and  /(-changing  RET  were 
accurately  modeled  with  four  parameters — one  for  dipole-dipole  collisions,  two  for  the  A J  scaling 
law,  and  one  for  the  cumulative  rate  of  /(-changing  collisions.  As  was  found  for  CH3F,  /-changing 
rotational  collision  rates  in  CH3C1  are  modeled  accurately  by  both  the  statistical  power  gap  (SPG) 
law  and  the  infinite  order  sudden  approximation  using  a  power  law  expression  for  the  basis  rates 
(IOS-P).  However,  in  contrast  to  CH3F,  where  all  IRMMDR  time  responses  for  /(-changing 
collisions  have  the  same  shape,  many  time  responses  of  CH3C1  states  populated  by  /(-changing 
collisions  contain  an  additional  early  time  feature  (ETF)  that  varies  with  pump  and  probe  states. 
Nonetheless,  a  simple  generalization  of  the  previously  reported  model  for  /(-changing  collisions  is 
shown  to  account  for  all  of  the  additional  features  observed  in  CH3C1.  Rather  than  observing  a  fixed 
temperature  for  /(-changing  collisions  as  was  the  case  for  CH3F,  the  temperature  is  found  to  be  a 
function  of  time  for  CH3C1.  Moreover,  the  two  new  parameters  this  adds  to  the  RET  model  are 
related  to  known  physical  quantities.  A  qualitative  argument  of  /(-changing  collisions  based  on  a 
classical  picture  is  offered  to  explain  the  difference  between  the  measured  /-  and  /(-changing 
state-to-state  rates  in  CH3C1. 


I.  INTRODUCTION 

Double  resonance  spectroscopy  using  infrared  pump  ra¬ 
diation  and  probe  radiation  in  the  millimeter-wave  to 
submillimeter-wave  range  (IRMMDR)  has  been  shown  to  be 
a  powerful  technique  for  the  study  of  both  rotational  and 
vibrational  state-changing  collisions  in  small  polyatomic 
molecules.1-7  In  our  earlier  time-resolved  IRMMDR  investi¬ 
gations  of  CH3F,  we  found  that  /-changing  collisions  (for 
u 3=  1  of  CH3F,  &J  =  n,  &K=0,  where  n  is  an  integer)  are 
modeled  well  by  both  the  statistical  power  gap  (SPG)  fitting 
law  and  the  infinite  order  sudden  approximation  assuming  a 
power  law  expression  for  the  basis  rates  (IOS-P).6  While 
dipole-dipole  rates  were  the  most  rapid,  quantum  number 
changes  of  up  to  |A/|  =  10  were  observed  in  the  data.  In 
contrast,  /(-changing  rotational  collisions  (for  u3  =  l  of 
CH3F,  A J  =  m,  &K=3n,  where  m  and  n  are  integers)  were 
not  found  to  be  selective  in  either  J  or  K  beyond  the 
symmetry  preservation  propensity  rule  discussed  by  Oka, 
A/(  =  3n  for  symmetric  top  molecules  of  C3l)  symmetry.8 
These  /(-changing  collisions  were  found  to  be  capable  of 
causing  large  changes  of  a  molecule’s  rotational  energy  in  a 
single  collision,  populating  rotational  states  in  proportion  to 
their  contributions  to  the  rotational  partition  function/  In 
other  words,  the  rotational  states  of  the  same  symmetry  type 
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and  vibrational  state  as  the  pumped  state  maintain  a  Boltz¬ 
mann  population  distribution  at  all  times.  We  refer  to  a  col¬ 
lisional  process  as  thermal  if  it  populates  a  group  of  states  in 
this  manner.  This  difference  between/-  and  /(-changing  col¬ 
lisions  was  observed  for  both  13CH3F  and  12CH3F.  A  com¬ 
puter  program  based  on  this  model  for  rotational  energy 
transfer  (RET)  shows  that  all  rotational  collisional  processes 
observed  in  CH3F  can  be  modeled  with  only  four  free 
parameters — one  for  dipole-dipole  collisions,  two  for  either 
the  SPG  or  IOS-P  law,  and  a  cumulative  rate  for  /(-changing 
collisions.6 

Using  time-resolved  IRMMDR,  we  have  studied  both 
rotational  and  vibrational  collision  processes  in  CH3C1. 
However,  the  present  paper  is  limited  to  the  discussion  of 
rotational  processes;  the  results  concerning  vibrational  pro¬ 
cesses  will  be  reported  later.  In  early  work  on  CH3C1  using  a 
similar  experimental  technique,  but  not  resolved  in  time  and 
using  a  limited  range  of  probe  frequencies,  Frenkel  et  aL 
found  that  |A/|  =  1  is  more  common  than  |A/|>1  in  A/(=0 
collisions,  but  that  in  A/(=- 3  collisions,  large  changes  in/ 
were  relatively  more  common.9  These  findings  are  consistent 
with  the  more  detail?-*  descriptions  of  RET  in  the  CH3F  stud¬ 
ies  of  Everitt  and  De  Lucia.5-6  The  purpose  of  the  present 
study  is  to  characterize  the  rotational  and  vibration*!  colli¬ 
sional  processes  of  CH3G  and  model  them  with  a  few  pa¬ 
rameters,  and  where  possible,  relate  these  parameters  to 
those  of  CH3F  and  fundamental  molecular  processes.  This 
objective  has  been  achieved  for  the  wide  range  of  states 
pumped  and  probed  in  this  experiment.  We  report  here  that 
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7 -changing  collisions  in  CHjCl  (for  o6=  1  of  CH3C1,  A7  =  /i, 
AAf  =  0,  A/=0,  where  n  is  an  integer  and  /  is  the  quantum 
number  for  angular  momentum  from  the  degenerate  u6  vi¬ 
brational  motion)  are  very  similar  to  those  in  CH3F.  How¬ 
ever,  a  difference  between  the  experimental  observations  of 
CH3C1  and  CH3F  led  to  a  generalization  of  the  model  for 
thermal  Af -changing  collisions  (for  u5=  1  of  CH3C1,  A J  =  m, 
A(Af- /)  =  3n,  where  m  and  n  are  integers  and  /  is  the  same 
as  above)  that  uses  a  small  number  of  parameters,  all  of 
which  are  physically  meaningful. 

Other  RET  studies  using  a  variety  of  experimental  tech¬ 
niques  have  been  performed  on  CH3F.  Oka  laid  out  the  basic 
propensity  rules  for  RET  in  symmetric  top  molecules  such  as 
CH3F  as  well  as  other  molecules  in  his  pioneering  micro- 
wave  double  resonance  experiments.8  In  ffee-jet  expansion 
experiments,  Douketis  and  co-workers  found  that  the  7  de¬ 
gree  of  rotational  freedom  of  CH3F  in  collisions  with  He 
cools  twice  as  effectively  as  that  of  Af.10  Schwendeman  and 
co-workers  have  studied  CH3F  in  a  series  of  frequency- 
resolved  infrared-infrared  double  resonance  experiments 
(IRIRDR)  observing  collisional  processes  in  addition  to  col¬ 
lisionless  effects.11'15 

Several  experimental  techniques  have  been  used  to  study 
RET  in  many  other  molecules  with  a  variety  of  collision 
partners.  A  large  body  of  experimental  and  theoretical  work 
with  diatomic  molecules  has  established  scaling  laws  that 
accurately  describe  RET  using  a  small  number  of 
parameters.16-23  The  rates  of  7-changing  RET  in  some  early 
studies  of  dipolar  molecules  [NH3,24  CDF^25  and  CH3F 
(Ref.  11)]  were  explained  using  only  A7  =  ±l  and  possibly 
AJ  =  ±2  collisions;  however,  recent  studies  have  shown  that 
large  changes  in  the  7  quantum  number  can  occur  in  single 
collisions  and  can  contribute  significantly  to  RET 
measurements.6'26'29  In  fact,  the  need  to  include  higher-order 
moments  (|A7|s=2)  in  NH3  was  suggested  as  early  as  1974. 30 
Orr  and  co-workers  found  7-changing  collisions  in  D2CO  are 
modeled  well  by  either  the  exponential  energy  gap  iaw 
(EGL)  or  the  SPG  law,  but  that  a  series  of  A7  =  ±l  and 
A7  =  ±2  collisions  alone  was  unable  to  reproduce  their 
measurements.26'31'32  Measurements  of  state-to-state  rates  in 
self-collisions  of  C2H2  have  shown  significant  rates  for  large 
changes  of  7. 33  Steinfeld  and  co-workers  reported  that 
7-changing  collisions  in  NH3  can  also  be  modeled  by  the 
EGL  fitting  law.29  Furthermore,  a  theoretical  treatment  by 
Parson34  35  showed  that  on  the  time  scale  of  collisions,  two 
spherical  top  molecules  studied  by  the  Steinfeld  group, 
13CD4  (Refs.  27  and  36)  and  SiH4  (Ref.  28)  can  be  treated  as 
symmetric  tops  as  a  result  of  temporary  rotational  distortion 
around  a  threefold  or  a  fourfold  symmetry  axis.  This  physical 
effect  results  in  “principal  channels”  for  RET  in  spherical 
tops  similar  to  7-changing  RET  in  symmetric  tops.  Hetzler 
and  Steinfeld  showed  that  RET  in  these  principal  channels 
can  also  be  modeled  by  the  EGL  fitting  law  with  changes  of 
7  up  to  |A7]  ~5  observed  for  SiH4.M  As  mentioned  above, 
Everitt  and  De  Lucia  showed  that  both  the  SPG  and  10S-P 
laws  accurately  model  7-changing  collisions  in  CH3F.6 

There  have  also  been  a  limited  number  of  experimental 
studies  of  Af -changing  rotational  collisions  in  nonlinear  mol¬ 
ecules,  ail  of  which  corroborate  the  strength  of  Oka’s  sym- 
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metry  preservation  rule  A K=mn,  where  m  is  the  order  of 
the  symmetry  axis  and  n  is  an  integer.8  The  study  of  Frenkel 
et  al.  gave  an  early  indication  that  Af -changing  collisions  in 
CH3C1  can  change  7  by  several  quantum  numbers.9  Several 
years  later  in  reports  based  on  time-resolved  measurements 
of  a  large  number  of  states,  De  Lucia  and  co-workers  re-lj 
ported  that  the  AAf  =  3n  process  in  CH3F  is  thermal  with  a 
cross  section  larger  than  the  gas  kinetic  collision  cross  - 
section.4'5  The  presence  of  a  AAf  =  3n  process  in  CH3F  has  '***'■• 
been  confirmed  by  Shin  et  al.  in  their  frequency  domain  IR-  "***’ ", 
IRDR  measurements.12  Harradine  et  al.  have  also  observed 
Af -changing  collisions  of  at  least  AAf=-15  in  self-collisicis“ 
of  CDFj.25  However,  their  results  do  not  indicate  that  the" 

Af -changing  process  in  this  oblate  top  is  thermal  in  both  J 
and  Af  since  they  found  that  only  a  narrow  distribution  of 
states  is  populated  by  Af -changing  collisions.  Orr  and  co-  ’ 
workers  have  also  discussed  Ka -changing  collisions  in  the  T 
asymmetric  top  D_C0.31  In  a  paper  concerned  primarily  with 
7 -changing  collision  rates,  they  modeled  population  loss  into 
states  of  different  Af  as  occurring  with  the  gas  kinetic  colli¬ 
sion  rate  for  only  AAC„  =  ±2  and  |A£|<10  cm-1,  while  all 
other  Af-changing  state-to-state  rates  were  zero.  Large 
changes  of  7  were  required  to  meet  the  limited-energy- 
change  requirement.  However,  in  more  recent  work,  they 
have  used  a  different  model  that  allows  large  changes  of  Af, 
AKa=±2n  to  conserve  symmetry,  but  only  zero  or  small 
changes  of  7  for  collisions  that  change  Af  ,32  Although  a  de¬ 
tailed  report  of  Af -changing  collisions  in  D2CO  has  not  been 
published,  it  appears  that  the  AKa  =  ±  2 n  process  is  not  ther¬ 
mal  in  this  near-prolate  top;  however,  large  changes  of  K  do 
occur.  Abel  et  al.  have  reported  state-to-state  rates  for 
Af -changing  collisions  in  NH3  self-collisions  and  in  collisions 
with  At.29  For  self-collisions,  a  comparison  of  7-specific 
=  — 3  collision  rates  with  AAf=0  collision  rates  shows  that 
the  Af -changing  collisions  are  less  probable  and  their  rates 
vary  more  slowly  in  7.  This  would  be  expected  since  they 
involve  higher-order  moments  of  the  intermolecular  poten¬ 
tial.  Abel  et  al.  obtained  reasonable  fits  to  the  individually 
measured  state-to-state  rate  constants  by  modeling  the  AAf 
=  —  3  collisions  using  the  energy  corrected  sudden  (ECS)  law 
and  basis  rates  that  vary  slowly  in  multipole  order.  Thus  for 
another  molecule.  Af -changing  collisions  have  been  shown  to 
be  less  state  selective  than  7-changing  collisions;  however, 
the  Af -changing  collisional  process  in  NH3  is  again  not  ther^ 
mal  as  in  CH3F  and  CH3C1. 


'-tC 


II.  BACKGROUND 

This  work  builds  on  an  understanding  of  the  collisional 
processes  of  CH3F.3'5'6,11'37  The  expectation  was  that  CH3C1 
would  exhibit  the  same  basic  collisional  processes  as  CH3F, 
although  perhaps  with  different  rates.  Since  the  CH3CI  model 
is  a  generalized  version  of  the  CH,F  model,  the  CH3F  work 
is  reviewed  in  the  following  section.  As  well  as  defining  the 
basic  components  of  the  model,  this  section  also  provides  a 
framework  for  comparing  the  rotational  processes  of  CHjF 
and  CH3CI. 


J.  Chem.  Phys.,  Vol.  100,  No.  8,  15  April  1994 


5666  ' 


Pape,  De  Lucia, and  Skatrud:  J  and  K  changing  in  CH3C1  collisions 


F  IG  1.  The  RET  model  used  for  CH,F  (Refs.  5  and  6).  Only  the  ground 
vibrational  slates  u,  =  l  and  u„  =  I  (the  lowest  energy  vibrational  modes) 
needed  to  be  included  to  accurately  model  the  IRMMDR  time  responses  for 
pump  intensities  that  were  used.  Examples  of  the  measured  collisional  pro¬ 
cesses  are  denoted  with  arrows.  While  only  two  types  of  ./-changing.  d/C =0 
processes  are  drawn,  they  were  all  modeled  using  a  scaling  law. 


A.  J -changing  collisions 

The  action  of  the  pump  in  this  double  resonance  experi¬ 
ment  is  to  populate  a  single  rotational  level  within  an  excited 
vibrational  state  with  molecules  from  a  single  rotational  state 
of  the  ground  vibrational  state.  Figure  1  shows  this  schemati¬ 
cally  for  1  CH,F.  In  the  earlier  Duke  Microwave  Laboratory 
work,1’6  two  coincidences  with  l2C1602  laser  lines  were 
used  extensively — 9P( 32)  pumps  13CH3F  population  from 
7=4.  K  =  3  of  the  ground  vibrational  state  to  ir3=  1 »  7=5, 
K  =  3.  and  9P( 20)  pumps  ‘“CHjF  population  from  7  =  12, 
K  =  2  of  the  ground  vibrational  state  to  o3=l,7=12,  K  =  2. 
As  the  pump  introduces  excess  population  in  the  excited  vi¬ 
brational  state,  several  collisional  processes  reestablish  rota¬ 
tional  thermal  equilibrium  within  1/3=1.  Likewise,  several 
vibrational  collisional  processes  operate  on  longer  time 
scales  to  equilibrate  population  among  the  vibrational  levels. 

The  fastest  of  the  collisional  mechanisms  is  the  dipoie- 
dipole  process.  The  well-known  selection  rules  for  this  pro¬ 
cess  are  A7=0.±l.  AK=0,  where  the  branching  ratios  are 
given  by  the  dipole  matrix  elements  for  a  symmetric  top38 


(7+ 1  )2- AT2 

*t-y+1-*dd  (_/  +  j )( 27  +  l )  ’ 
K 2 

and 

J2-Kz 

k ’  ~*dd 7(27+1  )  ‘ 


(1) 


where  kiA  is  the  dipole-dipole  collision  rate  and  the  matrix 
elements  have  been  summed  over  M .  In  addition,  a  detailed 
balance  factor  is  included  for  transitions  downward  in  en¬ 
ergy' 

Higher-order  collision  processes  that  change  the  7  state 
by  several  quanta  in  a  single  collision,  up  to  A7  =  - 10,  while 
leaving  K  unchanged  were  also  observed  in  CH3F.6  These 
collisions  are  characterized  by  selection  rules  A7  =  n, 
A/C=0,  where  n  is  determined  by  the  orders  of  the  interact¬ 
ing  multipole  moments.8  Examples  of  the  A7=2  process  are 
shown  in  Fig.  1.  Since  they  can  populate  states  several  7 
quantum  numbers  away  from  the  pumped  state  in  a  single 
collision,  the  effects  of  these  higher-order  collisions  are  easy 
to  observe  several  7  away  from  the  pumped  state  despite 
having  smaller  cross  sections  than  the  dipole-dipole  process; 
in  fact,  they  are  essential  for  reproducing  the  observed  popu¬ 
lation  flow  of  the  CH3F  system.  Multiple  dipole-dipole  col¬ 
lisions  alone  are  not  sufficient  to  account  for  the  observed 
IRMMDR  time  responses  of  states  several  7  away  from  the 
pumped  state.6  This  is  because  the  pseudorandom  walk  of 
the  dipole -dipole  process  in  this  one  dimensional  system 
requires  on  the  order  of  n~  collisions  to  move  population  n 
steps  in  7,  whereas  the  appropriate  higher-order  process  re¬ 
quires  only  one  collision  to  change  7  by  the  same  number. 

It  has  been  shown  that  both  the  infinite  order  sudden 
scaling  law  (IOS)  and  the  statistical  power  gap  fitting  law 
(SPG)  successfully  quantify  the  observed  7-changing  rate 
constants  in  CH3F.6  The  IOS  equation  for  A7  =  n,  AK=0 
collisions  in  CH3F  is16 

kjt->j  =[2J  f+  1  )expl — — — J  2  (2L  +  1) 


(7,  Jf 

\-K 


(2) 


where  (:::)  is  a  3-7  symbol  and  the  exponential  factor  ac¬ 
counts  for  detailed  balance,  7>  being  the  greater  of  7,  and 
Jf.  The  basis  rates,  which  could  be  derived  if  the  intenno- 
lecular  potential  were  known,  were  assumed  to  have  the 
power  law  form  kL^0  =  c[L(L  +  l)]~y  (thus  abbreviated 
lOS-P).  The  statistical  power  gap  law,  which  Smith  and  Prit¬ 
chard  showed  to  be  an  analytical  approximation  of  the  IOS-P 
scaling  law,22  was  calculated  using 


A£| 

\  kT  j 

kT 

(3) 


This  is  a  limiting  case  of  the  SPG  law  valid  in  the  regime 
where  |A7|/(7,+7^)  and  |A£|/kT  are  small  compared  to  1; 
both  are  true  for  the  measurements  of  CH3F. 


B.  K -changing  collisions 

Before  discussing  /C-changing  collisions,  we  would  like 
to  clarify  our  use  of  the  word  thermal.  In  this  paper,  thermal 
is  used  to  describe  both  state  population  distributions  and 
collisional  processes.  By  thermal  group  of  states,  we  mean 
that  the  population  within  these  states  can  be  described  by  a 
Boltzmann  distribution  with  a  single,  not  necessarily  ambi- 
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ent,  temperature.  Note  that  there  may  be  a  number  of  differ¬ 
ent  thermal  groups  and  that  the  total  population  in  a  group  is 
not  fixed  to  any  specific  value  and  may  even  vary  with  time. 
In  addition,  the  temperature  may  vary  with  time,  but  eventu¬ 
ally  both  the  temperature  and  the  total  population  of  the 
group  relax  to  ,heir  final,  ambient  values.  Often,  we  refer  to 
a  thermal  group  of  states  as  a  pool  of  states.  The  experimen¬ 
tally  determined  existence  of  thermal  pools  significantly  sim¬ 
plifies  the  modeling  of  population  flow  since,  for  each  ther¬ 
mal  group,  the  population  in  a  large  number  of  states  can  be 
described  by  just  two  variables.  For  example,  since  rotational 
processes  are  generally  faster  than  vibrational  processes,  the 
population  of  rotational  states  within  a  given  vibrational  state 
often  reaches  a  Boltzmann  distribution  (i.e..  becomes  ther¬ 
mal)  before  the  population  distribution  among  the  different 
vibrational  levels  reaches  equilibrium.  As  another  example 
for  a  symmetric  top,  the  rotational  states  of  the  same  vibra¬ 
tional  energy  and  symmetry  as  the  pumped  state  may  be¬ 
come  thermal  as  a  group  and  still  not  be  in  thermal  equilib¬ 
rium  with  the  states  of  the  other  symmetry  type,  even  though 
they  too  are  thermal  as  a  group. 

By  thermal  process,  we  mean  that  the  state-to-state  rates 
to  all  of  the  eligible  final  states  (states  satisfying  the  selection 
rules  of  the  process)  are  such  that  the  population  transferred 
into  these  final  states  is  distributed  according  to  a  Boltzmann 
distribution  with  a  single  temperature.  As  is  the  case  for  the 
temperature  of  the  thermal  pools,  the  temperature  for  the 
thermal  process  may  or  may  not  be  the  same  as  the  tempera¬ 
ture  of  the  cell.  Similar  to  the  thermal  pool  concept,  the 
motivation  for  modeling  processes  as  thermal  is  to  represent 
many  individual  state-to-state  rates  with  only  a  net  rate  and  a 
temperature.  Since  concepts  of  thermal  pools  and  rates  are 
justified  by  the  data,  we  use  them  to  greatly  simplify  the 
description  of  the  RET  processes,  thereby  yielding  a  rigor¬ 
ous,  yet  tractable,  model. 

In  addition  to  the  ./-changing  processes  discussed  above. 
A-changing  collisions  also  rapidly  depopulate  the  pumped 
state.  It  is  well  documented  that  because  of  nuclear  spin  sta¬ 
tistics.  changes  in  A  for  a  symmetric  top  molecule  with  C3u 
symmetry  are  restricted  to  multiples  of  three  on  the  time 
scales  of  experiments  such  as  this.8  Note  that  since  the  ±A 
levels  are  degenerate,  all  of  the  A#  3 n  states  are  coupled  by 
these  collisions  (e  g.,  K  =  1 — *  —  2).  Both  CH3F  and  CH3C1 
fall  into  this  symmetry  class.  Everitt  and  De  Lucia  showed 
the  AA=3n  process  in  CH3F  has  the  remarkable  character¬ 
istic  of  being  thermal.5  For  CH3F,  a  Boltzmann  distribution 
at  the  ambient  temperature  of  the  cell  (300  K)  adequately 
represented  the  distribution  of  the  AA  =  3 n  final  states.  This 
characteristic  was  tested  by  measuring  many  different  rota¬ 
tional  states  for  two  different  pump  coincidences  for  two 
different  isotopic  species.  UCH3F  and  13CH3F.  As  will  be 
discussed  later,  our  millimeter-submillimeter  wave  (mm/ 
submm)  probe  is  highly  sensitive  to  rotational  nonequilibria, 
so  the  experimental  evidence  is  compelling  that  the  states 
that  receive  population  through  the  A  =  3n  process  maintain 
a  thermal  distribution  of  population.  A-changing  RET  in 
CH,F  does  not  show  any  state  specificity  beyond  the  follow¬ 
ing  requirements:  the  vibrational  state  does  not  change;  the 
symmetry  type  does  not  change;  and  the  amount  of  popula¬ 


tion  a  state  receives  is  proportional  to  its  contribution  to  the 
rotational  partition  function.  This  last  requirement  is  in  sharp 
contrast  to  the  dipole-dipole  and  the  higher-order 
/-changing  collisionai  processes  which  are  distinguished  by 
/  selection  rules. 

The  thermal  A-changing  process  greatly  reduces  the 
complexity  of  the  description  of  the  CH3F  system.  Without 
such  a  process,  the  N  thermally  accessible  states  in  o3=l 
might  be  connected  by  up  to  N(N-  l)/2  separate  rates.  In¬ 
stead,  the  existence  of  the  thermal  AA=3n  process  allows 
all  of  the  rotational  states  of  different  A  than  the  pumped 
state  to  be  grouped  into  two  pools — the  pool  of  E  symmetry 
states  and  the  pool  of  A  symmetry  states. 

As  Everitt  pointed  out,37  it  is  not  clear  that  the  tempera¬ 
ture  of  the  distribution  of  AA=3n  final  states  should  neces¬ 
sarily  be  the  same  as  the  ambient  temperature  of  the  gas  as 
was  reported  for  CH3F.5  Based  on  an  energy  conservation 
argument,  it  might  be  expected  that  the  temperature  of  the 
distribution  of  final  states  should  be  different  from  the  cell 
temperature  since  the  rotational  energy  of  the  pumped  state 
is  in  general  different  than  EkT .  This  will  be  shown  to  be  the 
case  for  CH3C1. 

C.  Vibrational  state-changing  collisions 

The  previous  work  showed  that  for  CH3F,  some  vibra¬ 
tional  processes  operate  on  nearly  the  same  time  scale  as  the 
rotational  processes;  therefore,  they  cannot  be  ignored  when 
analyzing  experimental  RET  data.5  However,  since  the  focus 
of  this  paper  is  restrictea  to  the  rotational  processes,  we  will 
discuss  the  vibrational  processes  only  insofar  as  necessary  to 
deal  quantitatively  with  the  rotational  collisionai  mecha¬ 
nisms.  The  vibrational  processes  in  CH3C1  have  been  mea¬ 
sured  and  will  be  treated  further  in  a  future  paper. 

The  fastest  vibrational  process  in  CH3F  is  the  near¬ 
resonant  vibrational  swap  (F-swap)  mechanism 

CH3F(o3  =  1  )3  +  CH3F(o  =  0)s 

-CH3F(u  =  0)A  +  CH3F(t;3=l)8-A£,  (4) 

where  the  two  molecules  labeled  A  and  B  each  retain  their 
symmetry.  Here  the  vibrationallv  excited  molecule  A  ex¬ 
changes  a  unit  of  vibrational  energy  with  a  ground  state  mol¬ 
ecule  B.  If  the  colliding  molecules  are  of  opposite  symmetry, 
this  has  the  effect  of  transcending  the  symmetry  selection 
rule.  The  rate  for  this  process  in  CH,F  is  not  much  less  than 
the  rates  of  the  rotational  processes.5  7 

D.  Wall  collisions 

For  the  cell  size  and  gas  pressures  we  typically  use,  the 
most  important  mechanism  for  relaxing  the  fundamental  vi¬ 
brational  quantum  is  wall  collisions.  The  previous  CH3F 
analysis,  performed  using  a  2  cm  diameter  copper  tube, 
showed  that  a  wall  collision  has  a  high  probability  for  ther- 
malizing  a  molecule's  vibrational  energy.3  At  pressures  low 
enough  that  the  mean  free  path  becomes  long  compared  to 
the  radius  of  the  cell,  most  molecules  that  are  excited  by  the 
pump  move  to  the  wall  without  undergoing  any  velocity¬ 
changing  collisions.  (Note  that  we  can  still  observe  collision 
induced  IRMMDR  time  responses  at  these  low  pressures 
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IIG  2.  CH,C1  collision  rales  vs  pressure.  The  region  between  1  and  100 
mlbrr  is  the  pressure  range  for  the  measured  rotational  collision  data.  Note 
that  the  molecules'  motion  to  the  wall  changes  from  ballistic  to  diffusive  at 
a  pressure  of  about  5  mTorr  for  the  cell  of  4  82  cm  diameter 


since  many  of  the  state-to-state  rates  of  rotational  processes 
are  faster  than  the  gas  kinetic  collision  rate.)  Diffusion  be¬ 
gins  to  dominate  the  motion  to  the  wall  starting  at  a  pressure 
of  roughly  10  mTorT  for  room  temperature  CH,F  in  a 
2  cm  diameter  cell.  While  the  rate  of  diffusion  to  the  wail  is 
inversely  proportional  to  pressure,  the  rales  for  the  two-body 
collisional  processes  we  study  are  linear  in  pressure;  conse¬ 
quently.  it  is  easy  to  distinguish  relaxation  at  the  wall  from 
molecule -molecule  relaxation  processes  by  acquiring  data 
over  an  appropriate  range  of  pressures.  Figure  2  shows  the 
pressure  dependence  for  CH,C1  processes  in  a  4.82  cm  di¬ 
ameter  cell  used  in  the  present  experiment. 

III.  EXPERIMENT 

Although  the  experimental  setup  is  similar  to  that  re¬ 
ported  previously,4’*  the  apparatus  will  be  presented  here  to 
provide  a  description  of  recent  experimental  enhancements. 
Figure  3  shows  the  essential  elements  of  this  experiment.  A 
Q -switched  CO-  laser  coincidence  with  a  CH3C1  infrared 
transition  provides  a  nonthermal  population  distribution  in  a 
particular  rotational  state  of  the  oh=l  vibrational  state  of 


/ox  |  pzl  |  CO-  laser  iQ  switcnj — . 
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CHjCl.  A  mm/submm  spectrometer  simultaneously  probes  a 
single  rotational  transition  within  the  excited  vibrational 
state.  The  time-resolved  measurements  include  both  three- 
level  double  resonances,  where  the  pumped  state  is  also  one 
of  the  probed  states,  and  collisionallv  mediated  four-level 
measurements.  Using  a  model  of  me  collisional  processes 
and  the  measured  time-varving  absorption  lor  emission)  sig¬ 
nals.  a  nonlinear  least  squares  routine  then  fits  the  rates  of  the 
population  flow  processes. 

For  this  CH3C1  study,  a  higher  power  infrared  source 
was  used.  The  CO-  laser  is  a  50  W  Apollo  model  570  longi¬ 
tudinal  gas  flow  laser  that  we  modified  to  enable  Q  switching 
by  the  addition  of  an  intracavitv  CdTe  electro-optic  modula¬ 
tor.  This  system  allows  pulse  widths  to  be  adjusted  to  as 
narrow  as  0.5  /as  full  width  at  half-maximum  (FWHM).  Av¬ 
erage  powers  per  pulse  of  over  400  W  are  produced  for  the 
9PI26)  laser  line  at  a  1  kHz  repetition  rate.  The  diffraction 
grating  that  provides  single  line  oscillation,  while  allowing 
full  coverage  of  both  the  9  and  10  jum  bands  of  12C1602  is 
mounted  on  a  piezoelectric  translation  stage  to  allow  fine 
tuning  of  the  laser  cavity. 

The  mm.'submm  radiation  is  generated  using  harmonic 
generation  techniques  developed  in  the  Duke  Microwave 
Laboratory  by  King  and  Gordv  4  and  extended  to  the  1  THz 
range  by  Helminger.  Messer,  and  Dc  Lucia/11 41  A  phase  sta¬ 
bilized  35  GHz  klystron  is  impedance  matched  onto  a  point- 
contact  diode,  which  generates  radiation  at  integer  multiples 
of  the  klystron  frequency.  The  desired  harmonic  is  selected 
by  varying  the  size  of  the  multiplier's  output  waveguide  and 
by  tuning  the  impedance  matching  at  the  diode.  Interchang¬ 
ing  klystrons  and  output  frequency  cutoffs  allows  continuous 
coverage  of  the  mm/submm  spectrum  from  less  than 
100  GHz  to  greater  than  1  THz.  The  phase-lock  system  used 
in  this  experiment  yields  a  fractional  linewidth  of  better  than 
10”'  and  a  frequency  accuracy  of  better  than  10'8. 

The  body  of  the  diagnostic  cell  used  for  obtaining  the 
CH3C1  data  is  a  stainless  steel  pipe  1.0  m  long  by  4.82  cm 
diameter.  The  larger  diameter  of  this  new  cell  increases  the 
time  for  molecules  to  reach  the  cell  wall  and  thus  allows 
measurements  to  be  made  at  lower  pressures  where  colli¬ 
sional  rates  are  slower.  End  pieces  housing  dichroic  windows 
were  used  to  reflect  the  CO-  laser  radiation  and  to  transmit 
the  mrrusubmm  waves.  The  radii  of  both  the  pump  and  probe 
beams  are  expanded  to  about  80T  of  the  cross  sectional  ra¬ 
dius  of  the  cell.  After  the  beams  copropagate  the  full  length 
of  the  ceil  along  its  axis,  the  second  dichroic  window  splits 
the  beams;  the  laser  power  that  is  not  absorbed  by  the  gas  is 
monitored  at  a  Scientech  362  power  meter,  while  the  mm/ 
submm  radiation  is  directed  to  the  spectrometer's  broadband 
detector.  The  measurements  reported  here  were  made  at  pres¬ 
sures  of  I- 1000  mTorr  and  at  room  temperature. 

A  significant  improvement  to  the  experiment  over  that 
used  in  previous  papers  from  this  laboratory  is  the  ability  to 
distribute  both  pump  and  probe  radiation  smoothly  over  the 
cross  section  of  the  diagnostic  cell.  This  is  achieved  by  using 
large  aperture  dichroic  windows  to  combine  the  radiation 
rather  than  the  former  scheme  where  the  laser  beam  was 
focused  through  a  pinhole  in  the  cell  wall.  The  dichroic  win¬ 
dows  were  fabricated  of  a  multilayer  coating  of  Ge  and  ZnS 
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on  a  low-water  fused  silica  substrate.  At  the  45°  incidence 
used  in  the  experiment,  the  windows  have  a  reflectivity  of 
greater  than  92%  across  the  9  /j.m  band  of  the  CO:  laser  and 
94.3%  at  the  9.6  /xm  9P(26)  line.  The  laser  power  that  is  not 
reflected  at  the  dichroic's  coated  surface  is  absorbed  by  the 
substrate;  the  spectrometer's  broadband  detector  did  not 
measure  spurious  laser  radiation.  We  measured  these  win¬ 
dows  to  be  no  worse  than  75%  transmissive  to  both  150  and 
300  GHz  radiation  at  45°  incidence. 

The  mm/submm  detector  is  an  InSb  hot  electron  bolom¬ 
eter  operated  at  the  He-cooled  temperature  of  about  1.7  K. 
The  detector  bandwidth  is  about  1  MHz.  Improved  signal 
averaging  is  now  provided  by  two  computer-aided  measure¬ 
ment  and  control  ICAMAC)  modules,  the  dsp  Technologies 
2030S  transient  digitizer,  and  4101  averaging  memory.  This 
system  captures  every  waveform  at  the  typical  repetition  rate 
of  1  kHz.  In  most  cases,  60  000  waveforms  were  summed  for 
a  total  integration  time  of  1  min.  Most  of  the  signals  were 
averaged  over  a  100  fis  window  at  a  digitization  rate  of 
0.133  /us:  however,  a  300  [is  window  and  0.333  /xs  rate  were 
used  for  the  slowest  IRMMDR  time  responses. 

Time-resolved  IRMMDR  is  particularly  well-suited  for 
studying  rovibrational  collisional  kinetics  because  of  the  sen¬ 
sitivity  of  the  mm/submm  rotational  absorption  (and  emis¬ 
sion)  measurements,  particularly  to  nonthermal  rotational 
population  distributions.  While  the  experiment  is  sufficiently 
sensitive  to  measure  time  responses  for  states  populated  ther¬ 
mally  (with  a  signal  to  noise  ratio  of  better  than  10:1  for  1 
MHz  bandwidth  and  60  s  integration  for  many  pump- 
symmetrv  transitions),  the  experiment  is  much  more  sensi¬ 
tive  to  states  populated  nonthermally.  The  large  nonthermal 
signals  are  in  some  cases  due  to  the  large  pump-induced 
deviations  of  the  rotational  populations  from  thermal  distri¬ 
butions.  In  addition,  rotational  absorption  probing  of  this  re¬ 
gion  of  the  spectrum  is  intrinsically  very  sensitive  to  even 
small  deviations  from  thermal  population  distributions.  This 
sensitivity  is  crucial  to  the  verification  of  the  thermal  nature 
of  some  of  the  collisional  processes.  The  sensitivity  to  non¬ 
thermal  deviations  is  illustrated  by  the  following  example;  at 
212  GHz,  there  is  only  a  3.4%  thermal  population  difference 
between  probed  rotational  states,  and  this  difference  pro¬ 
duces  the  easily  observed  absorption  of  0.5%  in  the  excited 
vibrational  state  (even  though  rotational  states  in  the  oft=l 
vibrational  state  have  only  1/140  as  much  population  as 
those  in  the  ground  vibrational  state).  Therefore,  if  the  popu¬ 
lation  of  one  state  deviates  from  thermal  by  a  factor  of  only 
0.1%.  then  the  absorption  will  change  by  0.015%,  which 
results  in  a  signal  to  noise  of  about  5:1  for  our  typical  time- 
resolved  bandwidth  of  1  MHz  and  60  s  integration. 

IV.  SPECTROSCOPIC  CONSIDERATIONS 

Table  1  contains  the  pump  and  probe  transitions  used  in 
this  work.  The  IRMMDR  time  responses  used  for  fitting  the 
collision  rates  were  taken  when  the  laser  was  tuned  to  maxi¬ 
mum  power  and  the  mm/submm  probe  was  tuned  to  the  line 
center :  however,  many  other  time  responses  were  acquired 
with  the  pump,  probe,  or  both  offset  in  frequency.  The  iso¬ 
topic  species  of  CH,C1  were  present  in  natural  abundances 
76%  35C1  and  24%  37C1.  For  the  large  part,  we  both  pumped 


and  probed  CH,!<iCl:  however,  we  did  spot  check  all  three 
other  permutations — pump  CHi’5CI.  probe  CH,37C1;  pump 
CH337CI,  probe  CH,VSCI:  and  pump  CH,i7Cl.  probe  CH337C1. 
A  description  of  the  isotopic  species  study  will  be  reported  in 
a  future  paper.  In  contrast  to  those  used  with  CH?F,  none  of 
the  CO;-CH,C)  coincidences  used  in  this  study  displayed 
any  behavior  indicating  multiple  coincidences  for  a  given 
pump  line.  For  '*CH3F,  the  9P(20)  pump  line  put  small 
populations  in  AT  =  1  and  £  =  ()  of  7  =  11  in  addition  to  the 
main  pumped  state  7=11.  K  =  2. 

Several  spectroscopic  circumstances  conspire  to  make 
this  CH,C1  study  more  challenging  than  that  of  the  CH3F 
work,  two  of  these  being  methyl  chloride's  hyperfine  struc¬ 
ture  and  /-type  doubling  in  its  pumped  vibrational  state.  Fig¬ 
ure  4  shows  the  rotational  energy  levels  of  vh=  1  for  CH3CI. 
Chlorine's  large  quadrupole  moment  causes  hyperfine  split¬ 
ting  that  is  resolved  over  most  of  our  working  pressure 
range,  making  the  CH,C1  spectrum  dense.  In  addition,  /-type 
doubling  is  an  issue  in  the  CH,C1  work  because  the  C02 
laser  coincidences  with  CH^Cl  involve  transitions  from  the 
ground  state  to  the  doubly  degenerate  t/h=l  bending 
state.  K  is  the  total  angular  momentum  about  the  molecular 
symmetry  axis  consisting  of  pure  rotational  angular  momen¬ 
tum  K'  and  a  contribution  from  the  vibrational  motion  /.  The 
overall  symmetry  of  a  rotational  state  is  determined  by 
K'=K  —  l,  where  A  symmetry  slates  consist  of  those  with  K' 
equal  to  a  multiple  of  three  and  the  £  symmetry  states  con¬ 
sist  of  the  remaining  states.  The  svmmetrv  preservation  pro¬ 
pensity  rule  that  applies  to  RET  in  the  /-type  doubled  u6=l 
vibrational  manifold  of  CH3C1  is  MK~  l)  =  3n.  where  n  is 
an  integer. 

Both  hyperfine  structure  and  /-type  doubling  increases 
the  density  of  the  CH3C1  rotational  spectrum  compared  to  the 
spectrum  for  o3=l  in  CH,F.  The  hyperfine  structure  also 
reduces  the  absorption  coefficients  of  the  resolved  lines.  The 
weaker  pump  and  probe  absorption  coefficients  for  t/6=l  of 
CH3C1  made  necessary  many  of  the  experimental  improve¬ 
ments  mentioned  above.  Although  hyperfine  structure  and 
/-type  doubling  contribute  to  the  complexity  of  the  CH3C1 
spectrum,  they  also  provided  additional  rotational  states  to 
probe  and.  therefore,  the  possibility  of  gaining  additional  in¬ 
formation  about  the  rotational  collision  processes. 

Another  significant  difference  with  CH3F  is  that  CH3C1 
has  a  convenient  C02  pump  coincidence  that  populates  a 
state  with  rotational  energy  greater  than  room  temperature 
thermal  energy.  The  CH,F  pump  lines  used  in  the  earlier 
study  correspond  to  energies  of  64  cm  1  for  the  9P(32)  co¬ 
incidence  with  i3CH3F  and  to  148  cm  1  for  9P(20)  with 
i2CH3F,  both  of  which  are  below  the  205  cm"1  room  tem¬ 
perature  thermal  energy.  In  contrast,  the  CH3C1  pump  coin¬ 
cidences  are  85  cm  for  9P(26)  and  285  cm  1  for  9/?(12). 
As  will  be  discussed  below,  the  energy  difference  between 
the  pumped  states  in  CH3C1  aided  in  the  analysis  of  its  col¬ 
lision  processes. 

V.  EXPERIMENTAL  RESULTS 

A.  Explanation  of  IRMMDR  time  response  data 

Based  on  experimental  measurements,  the  IRMMDR 
time  responses  of  the  transitions  we  measured  may  be 
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TABLE  I.  Pump  and  probe  transitions  used  in  the  iRMMDR  measurements  of  CH>"CI. 


Pump  coincidences 


Frequency 

CO;  line  (cm1) 


Ground  slate  uh=  L 

J  K  J  K  I 


Frequency  offset 
(MHz) 


9/1(12)  1073.278  11  6  12  7  1  -30 

9P(  26)  1041.280  A  3  6  4  1  20 

Probe  Transitions 


Transition  Frequency' 


Vib 

J-J' 

K 

/ 

Symmetry 

(GHz) 

v6=l 

3-4 

-  I 

-1 

A 

105.8915 

‘'#  =  1 

3-4 

1 

1 

A 

106.0079 

u6=l 

3-4 

2 

-1 

A 

105.9353 

u6=l 

4-5 

- 1 

-1 

A 

132.3594 

v„=\ 

4-5 

i 

1 

A 

132.5051 

t>6=l 

4-5 

2 

-1 

A 

132.4103 

u6=l 

4-5 

4 

1 

A 

132.4213 

V6=l 

5-6 

- 1 

-1 

A 

158.8268 

U6=l 

5-6 

i 

1 

A 

159.0016 

U6=l 

5-6 

-1 

A 

158.8910 

Vft=l 

5-6 

4 

1 

A 

158.8977 

”6=1 

5-6 

5 

-1 

A 

158.8130 

U6=l 

6-7 

-1 

-1 

A 

185.2913 

t>s=l 

6-7 

1 

1 

A 

185.4952 

V„=l 

6-7 

? 

-1 

A 

185.3643 

U6=l 

6-7 

4 

1 

A 

185.3725 

V6=l 

6-7 

5 

-1 

A 

185.2888 

U#=l 

7-8 

-1 

-1 

A 

211.7528 

t>6=l 

7-8 

1 

1 

A 

211.9857 

U6=l 

7-8 

1 

-1 

A 

211.8366 

U6=l 

7-8 

4 

1 

A 

211.8448 

u6=l 

7-8 

5 

-1 

A 

211.7428 

Ut=l 

7-8 

■7 

I 

A 

211.7660 

u#=l 

8-9 

-  1 

-1 

A 

238.2109 

U6=l 

8-9 

1 

1 

A 

238.4728 

U6=l 

8-9 

-1 

A 

238.3055 

U(,=  l 

8-9 

1 

1 

E 

238.3399 

U6=l 

8-9 

3 

-1 

E 

238.2788 

U4=l 

8-9 

4 

-1 

E 

238.2447 

^6=1 

8  9 

4 

1 

A 

238.3137 

U6=l 

8-9 

5 

-1 

A 

238.2018 

U6=l 

8-9 

6 

-1 

E 

238.1553 

Vft—  1 

8-9 

7 

1 

A 

238.2213 

u6=l 

9-10 

-1 

A 

264.7700 

u6=l 

9-10 

4 

1 

A 

264.7787 

u6=l 

10-11 

2 

-1 

A 

291.2302 

^6=1 

10-11 

4 

1 

A 

291.2396 

U,=  l 

4-5 

3 

A 

131.7527 

u3=l 

8-9 

0 

A 

237.1448 

Uj=l 

8-9 

1 

E 

237.1412 

u,  =  l 

8-9 

2 

E 

237.1305 

u,=  l 

8-9 

3 

A 

237.1124 

u,  =  l 

8-9 

4 

E 

237.0876 

'Here  hyperfine  components  are  resolved  at  20  mTorr.  the  rotational  frequencies  were  chosen  to  minimize 
spurious  signals  from  nearby  transitions. 


grouped  into  two  categories — those  where  the  states  in¬ 
volved  in  the  transition  have  a  nonthermal  population  distri¬ 
bution  and  those  involving  states  that  maintain  a  thermal 
population  distribution. 

Nonthermal.  For  the  first  several  microseconds  follow- 
“g  the  pump  pulse,  the  transitions  involving  states  of  the 
*®me  K  and  l  as  the  pumped  state  displayed  strong  nonequi¬ 


libria  due  to  the  fast  dipole-dipole  and  other  significant 
multipole-multipole  collisional  processes.  (We  refer  to  a  set 
of  J  states  with  the  same  K  and  l  and  within  the  same  vibra¬ 
tional  state  as  a  “AT//  stack.”)  The  signal  to  noise  estimation 
in  the  Experiment  section  above  shows  the  high  sensitivity  to 
these  nonthermal  IRMMDR  time  responses  in  the  pumped 
K/l  stack. 
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3  4 

K'=K-l 


FIG.  4.  Rotational  energy  levels  of  the  CHj]5C1  ufi  =  l  vibrational  manifold,  which  experiences  /-type  doubling.  Hyperfine  components  of  a  state  have  been 
reduced  to  a  single  value  by  talcing  a  weighted  average.  Kll  stacks  are  paired  by  K '  with  the  /  =  -  1  stack  on  the  left  and  the  /  =  1  stack  on  the  right.  The  two 
states  that  are  populated  bv  the  CO-  pump  are  labeled  9P(26)  and  9/?(12).  Transitions  for  which  we  have  acquired  1RMMDR  time  responses  are  marked  with 
vertical  lines  The  A  symmetry  states,  the  pump  symmetry  for  both  coincidences,  are  shown  wider  than  the  £  symmetry  states. 


Thermal.  Time  responses  involving  all  other  u6=  1  states 
belonged  to  one  of  two  thermal  pools — the  A  symmetry  pool 
and  the  E  symmetry  pool.  As  stated  before,  all  states  within 
a  pool  remain  in  thermal  equilibrium  among  themselves, 
even  though  both  the  temperature  and  total  population  of  the 
pool  are  free  to  change  in  time.  However,  if  the  pool  tem¬ 
perature  remains  fixed  at  the  temperature  of  the  cell,  then  all 
of  the  time  responses  of  a  pool  have  the  same  amplitude- 
normalized  shape  as  was  observed  in  CH3F.5 

The  probe  monitors  a  difference  in  population  between 
two  adjacent  J  states.  We  checked  for  M  selectivity  in  the 
observed  time  responses  involving  the  pumped  state  by  mak¬ 
ing  measurements  with  the  linear  pump  and  probe  polariza¬ 
tions  parallel  and  perpendicular  to  each  other;  however,  we 
did  not  observe  any  systematic  difference  in  the  time  re¬ 
sponses  as  has  been  observed  in  H2C0.42  Consequently,  the 
model  assumes  all  M  levels  are  populated  equally.  For  our 
IRMMDR  data,  the  zero  absorption  level  of  a  time  response 
corresponds  to  an  ambient  equilibrium  population  difference 
between  the  two  states;  consequently,  absorption  measure¬ 
ments  are  relative  to  the  signal  from  'he  thermal  “back¬ 
ground”  population  since  the  pump  does  not  significantly 
heat  the  gas.  If  the  population  difference  between  the  lower 
and  upper  states  decreases,  the  delector  measures  a  negative 
signal  which  corresponds  to  the  decreased  absorption  (more 
probe  power  arrives  at  the  detector).  However,  since  we  are 
concerned  with  the  population  differences  that  deviate  from 


0  5  10  15  20 

time  (ms) 


FIG.  5.  IRMMDR  time  responses  of  several  transitions  in  the  pumped  Kll 
stack.  The  state  populated  by  the  9P(26)  line  is  J  - 6,  K= 4,  /  =  1  and  the 
state  populated  by  the  9/?(12)  line  is./ =  12.  AC  =  7,  /  =  1.  The  data  (dots)  were 
taken  at  20  mTorT.  The  (its  for  both  pump  coincidences  (solid  lines)  were 
performed  simultaneously  for  these  time  responses  and  44  others  at  different 
pressures  using  the  IOS-P  scaling  law  for  all  A/f =0  collisions. 
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|9P(26)| 


|9R(12)| 


J  =  8-9 


J=  7-8 


J  =  6-7 


4-5 


K*  1,  /  =  1 


K»2,/  =  -1  K»4,r«1 


K=5,  r  =  -1 


K*  7.  /  *  1 


FIG  6.  Measured  CH,G  u«=  t,  4  symmetry  IRMMDR  time  responses  for  several  rotational  states  at  30  mTorr  with  fit  time  responses  for  two  different  pomp 
lines.  9PI26)  and  9R(12).  Each  vertical  axis  is  the  normalized  millimeter-wave  absorption  and  each  horizontal  lime  axis  covers  60  /js.  The  graphs  are 
arranged  in  the  same  order  as  the  A  symmetry  energy  levels  shown  in  Fig.  4.  The  dots  are  the  data,  the  solid  lines  are  the  fit  time  responses,  and  the  dashed 
lines  indicate  the  equilibrium  absorption  levels.  For  9/>(26),  the  pumped  KH  stack  is  K =4,  /  =  l  and  for  9R(12),  it  is  K  =  7.  f=l.  The  fitted  time  responses 
for  7 -changing  collisions  were  modeled  using  the  IOS-P  model. 


the  equilibrium  distribution,  we  refer  to  any  decreased  ab¬ 
sorption  as  emission,  even  though  the  net  effect  may  just  be 
reduced  absorption.  Furthermore,  since  the  thermal  back¬ 
ground  population  is  approximately  constant  in  time,  only 


the  pumped  molecules  undergo  state-changing  collisions  in 
the  model.  Accordingly,  the  thermal  background  population 
is  not  included  in  the  definitions  of  the  pools  for  modeling 
purposes,  and  when  we  refer  to  the  Boltzmann  distribution  of 
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population  in  a  pool,  it  is  only  the  pumped  molecules  that 
have  not  yet  been  vibrationally  the  malized  that  we  consider. 

B.  A  comparison  with  CH3F 

As  expected,  there  were  many  similarities  between  the 
experimental  results  for  CH3F  and  CH3C1.  There  were  sev¬ 
eral  types  of  IRMMDR  time  responses  in  both  molecules — 
nonthermal  time  responses  for  J -changing  collisions,  and  the 
time  responses  of  several  thermal  pools.  However,  the  time 
responses  for  CH3C1  states  populated  by  K-changing  colli¬ 
sions  showed  some  significant  differences  from  the  CH3F 
data. 

Figure  5  shows  the  IRMMDR  time  responses  of  several 
transitions  in  the  pumped  K/l  stacks  of  CH3C1.  The  behavior 
of  these  states  is  very  similar  to  those  of  the  pumped  K  stack 
in  CH3F.  The  largest  and  fastest  IRMMDR  time  responses 
come  from  the  three-level  double  resonance  measurements. 
Time  responses  in  the  same  K/l  stack  as  the  pumped  state 
become  increasingly  smaller  and  slower  with  increasing 
\\J\.  For  the  9/? (12)  pump,  all  of  the  measured  transitions  in 
the  pumped  K/l  stack  were  below  the  J=\2  pumped  state; 
consequently,  all  exhibited  emission. 

The  IRMMDR  time  responses  for  CH3C1  states  popu¬ 
lated  by  vibrational  processes  are  also  similar  to  those  of 
CH?F;  however,  the  F-swap  process  in  CH3C1  is  slower 
than  the  u3,  Inswap  process  in  CH3F.  This  is  due  to  the 
smaller  transition  moment  for  bending  mode  transitions  in 
u6=l  of  CH3C1. 43-44  Additionally,  the  u3=l  states  are  popu¬ 
lated  by  the  nonresonant  p3«— >v6  vibrational  process,  and  the 
V—T/R  process  is  an  important  vibrational  relaxation  pro¬ 
cess  at  higher  pressures  where  the  diffusion  rate  to  the  wall  is 
slower. 

While  the  IRMMDR  time  responses  of  the  above  CH3C1 
transitions  are  very  similar  to  those  of  CH3F,  the  time  re¬ 
sponses  of  the  states  populated  by  K-changing  collisions 
li  e.,  the  u6=l,  A  states)  are  different.  Figure  6  shows  a 
number  of  CH3C1  time  responses  for  the  t/6=l,  A  states  us¬ 
ing  both  pump  lines.  Similar  data  were  acquired  for  other 
pressures  from  1  to  100  mTorr.  The  following  are  some  of 
ihe  salient  features  of  these  time  responses:  (1)  K-changing 
collisions  are  fast.  As  was  the  case  in  CH3F,  K-changing 
collisions  occur  at  a  rate  comparable  to  gas  kinetic  collisions. 
Furthermore,  the  cumulative  rate  for  K-changing  collisions  is 
within  an  order  of  magnitude  of  the  rate  of  the  dipole-dipole 
process.  i2)  IRMMDR  time  responses  are  different  for  the 
9P(26)  and  9R(12)  pump  lines.  As  seen  in  9R(\2)  time  re¬ 
sponses  of  Fig.  6,  the  measurements  with  the  9£(12)  pump 
coincidence  are  similar  to  those  of  CH3F,  in  that  all  of  the 
time  responses  have  the  same  shape  to  within  the  noise  level. 
However,  the  measurements  made  using  the  9P(26)  pump 
line  reveal  an  additional  early  time  feature  (ETF)  that  was 
not  observed  in  the  CH3F  data.  An  ETF  is  the  fast  rise  and 
fall  seen  in  a  9P(26)  time  response  at  early  time,  especially 
at  lower  K',  that  is  absent  in  the  same  transition  when  using 
the  9R(12)  pump  line.  (3)  IRMMDR  time  responses  are  ther¬ 
mal  in  J.  This  is  true  for  the  data  of  both  pump  coincidences. 
The  sharp  nonthermal  features  found  in  the  pumped  K/l 
stack  data  of  Fig.  6  are  not  present  in  the  time  responses  of 
the  transitions  that  are  populated  by  K-changing  collisions. 
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FIG.  7.  IRMMDR  time  responses  of  CH,G  A  -symmetry  states  in  the 
vibrational  manifold  for  J  =  3-4  to  8-9.  The  time  responses  are  offset  ver¬ 
tically  for  legibility.  All  data  ait  from  f(=l.  /=1  using  the  9P(26)  pump 
line  and  a  gas  pressure  of  50  mTorr.  The  spike  at  early  time  in  J  =6-1  is  due 
to  the  a.c.  Stark  effect. 


Even  though  the  time  responses  of  9£(26)  display  a  distinc¬ 
tive  ETF,  those  within  a  nonpumped  K/l  stack  have  the  same 
shape.  Figure  7  further  demonstrates  this  point.  Given  the 
mm/submm  probe’s  sensitivity  to  nonthermal  population  dis 
tributions  between  adjacent  J  states,  these  data  show  that 
K-changing  RET  for  a  given  K/l  stack  is  very  close  to  ther 
mai.  (4)  For  9P(26),  the  ETF  is  larger  in  low  energy  states. 
At  early  time,  more  pump  population  that  has  undetgone 
K-changing  collisions  is  found  in  the  K'  =  0  states  (K=-l, 
/=  —  1  and  K  =  1,  1=1)  than  is  found  in  K'=  6  states  (K=5, 
/=- 1  and  K  =  7,  1=1).  Since  the  changes  in  quantum  num¬ 
ber  from  the  pumped  to  probed  state  are  similar  in  both  cases 
and  the  energies  of  the  K'=6,  K/l  stacks  are  significantly 
higher  than  those  of  K'=0  and  K'=  3  (see  Fig.  4),  the  ETF 
apparently  has  a  stronger  dependence  on  the  state’s  rotational 
energy  than  on  the  state’s  angular  momentum.  (5)  For 
9R(12),  a  wide  range  of  states  are  populated  thermally.  Fig 
ure  6  shows  that  K-changing  RET  in  9/?(  12)  yields  the  same 
IRMMDR  time  responses  for  many  different  states.  Picking 
the  extremes  of  our  measurements,  the  J=4,  K=l,  /  =  1  state 
(rotational  energy  11  cm’1)  is  populated  at  about  the  same 
rate  as  J  =9,  K=5,  l  =  -\  (rotational  energy  172  cm  ').  This 
is  a  clear  demonstration  that  K-changing  RET  is  not  state 
specific.  It  transfers  population  over  a  change  of  AJ=8 
AK=6,  A£  =  274  cm’1  at  a  rate  comparable  to  that  of  A J 
=  3,  AK  =  2,  A£  =  113  cm’1. 


VI.  ANALYSIS 

A.  Generalized  model  for  RET 

The  initial  CH3CI  analysis  was  performed  using  an 
adapted  version  of  the  computer  program  developed  for  the 
CHjF  analysis.  The  CH,F  model  is  represented  schemati¬ 
cally  in  Fig.  1.  This  successfully  modeled  collisionalwmgy 
transfer  in  CH3C1,  except  that  the  state-dependent  ETF  of  the 
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9P(26)  IRMMDR  time  responses  could  not  be  reproduced 
since  the  CH3F  model  requires  that  all  time  responses  be  the 
same  shape  for  all  A  states.  However,  a  straightforward  gen¬ 
eralization  of  the  AT-changing  RET  process  as  modeled  in 
CH3F  fully  describes  all  of  the  features  observed  in  the 
CHjCI  data,  including  the  ETF.  In  the  CH3F  model. 

-changing  collisions  distribute  the  population  among  the 
various  7  states  according  to  a  Boltzmann  distribution  at  a 
fixed  temperature,  specifically,  the  temperature  of  the  cell. 
For  CH3C1,  while  AT -changing  RET  still  populates  final  states 
according  to  thermal  Boltzmann  statistics,  the  temperature  of 
the  population  transferred  by  these  collisions  is  initially 
some  nonambient  temperature  that  subsequently  relaxes  to 
the  equilibrium  temperature.  The  time  dependence  of  this 
temperature  is  assumed  to  be 

r(r)  =  rf  +  (T0-7>-*™,  (5) 

where  Te  is  the  equilibrium  ambient  temperature  of  the  gas, 
Ta  is  the  initial  temperature  of  the  pool,  and  (fcTR  is  the  rate 
constant  of  the  temperature  relaxation. 

A  simple  argument  based  on  energy  conservation  moti¬ 
vates  this  modification  of  the  model.  Figure  4  shows  the 
energies  of  the  pumped  states  relative  to  room  temperature 
thermal  energy  EkT.  At  equilibrium,  both  partners  in  a  colli¬ 
sion  come  from  a  thermal  distribution  of  rotational  and  trans¬ 
lational  energies.  When  the  pump  pulse  populates  a  specific 
rotational  state  in  the  u6=  1  vibrational  state,  most  of  these 
pumped  molecules  still  collide  with  molecules  that  are  ther¬ 
mal  since  fewer  than  1 %  of  all  molecules  are  excited  by  the 
pump.  However,  the  pumped  molecules  are  not  thermal  in 
their  rotational  energy.  (They  may  also  have  a  nonthermal 
translation  energy  which  will  be  discussed  later.)  When  the 
first  AT-changing  collision  occurs  for  a  molecule  pumped  by 
the  9P(26)  pump  line,  most  of  the  nonthermal  population  is 
in  only  a  few  J  states  of  AT  = 4,  /  =  1  since  the  dipole-dipole 
process  is  limited  in  range  of  J  and  is  only  marginally  faster 
than  the  cumulative  rate  of  AT-changing  collisions.  The  rota¬ 
tional  energies  of  7=5-7  are  between  79  and  91  cm-1. 
While  its  collision  partner  is  almost  always  a  thermal  ground 
state  molecule,  the  initial  state  of  the  probed  molecule  is 
prepared  by  the  pump  to  have  one  of  these  low  rotational 
energies.  Thus  the  average  rotational  energy  available  in  the 
molecule’s  first  AT -changing  collision  is  less  than  Ekr,  so  it  is 
not  surprising  that  the  rotational  distribution  of  final  states 
from  these  collisions  is  at  predominately  lower  energies. 

For  the  9P(26)  pump  line,  the  result  of  subsequent 
AT-changing  collisions  is  to  raise  the  temperature  of  the  A 
symmetry  pool,  so  that  the  final  temperature  of  the  pool  is 
the  same  as  the  cell  temperature.  As  we  argued  above,  the 
first  AT -changing  collision  of  each  pumped  molecule  leaves  it 
in  a  lower  energy  state  on  average.  However,  the  data  indi¬ 
cate  that  the  AT-changing  RET  process  in  CH3C1  is  not  state 
selective  beyond  the  requirement  that  the  initial  temperature 
of  the  distribution  be  lower  than  the  cell  temperature;  the 
final  state  of  an  individual  rotational  collision  can  be  several 
7  and  AT  away  from  the  initial  state  and  can  have  a  very 
different  energy.  Since  almost  ail  AT -changing  collisions  in¬ 
volve  thermal  ground  state  molecules,  AT -changing  collisions 
are  efficient  at  uniformly  raising  the  temperature  of  the  cold 
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FIG.  8.  Several  J  =8-9  IRMMDR  time  responses  of  K -changing  collisions 
for  9/>C26)  and  9fl(  12)  at  50  mTorr.  One  set  of  fitted  time  responses  (dasbed 
lines)  shows  the  effect  of  fixing  the  temperature  of  the  /(-changing  RET 
process  to  the  cel)  temperature,  while  the  other  set  (solid  lines)  comes  from 
the  generalized  model  where  the  temperature  of  transferred  population  de¬ 
cays  to  the  cell  temperature  after  starting  at  a  nonambient  initial  tempera¬ 
ture.  The  discrepancy  in  the  first  microsecond  of  the  9/>(26)  K= 2,  /=-l 
time  response  is  due  to  the  Lorentzian  tail  of  the  much  stronger  signal  at 
K= 4,  /=  1. 


A  state  distribution.  In  other  words,  the  temperature  of  the  A 
pool  relaxes  towards  the  cell  temperature  at  a  rate  slightly 
slower  than  the  cumulative  rate  of  AT-changing  collisions  be¬ 
cause  it  takes  a  few  of  these  collisions  to  cause  the  relax¬ 
ation. 

While  the  above  argument  was  made  for  the  9P(26) 
pump  line,  the  fact  that  the  ETF  is  not  discernible  in  the 
9/?(12)  data  is  also  consistent  with  this  temperature-relaxing 
picture.  Whereas  AT -changing  RET  for  the  9P(26)  pump  ini¬ 
tially  concentrates  the  A  pool  population  in  the  lowest  rota¬ 
tional  levels,  AT-changing  collisions  for  the  9/?(12)  pump  ini¬ 
tially  spread  the  pumped  population  over  a  larger  number  of 
A  states  due  to  the  higher  rotational  energy  of  the  state 
pumped  by  the  9A?(12)  coincidence.  Just  as  for  the  case  for 
the  9P(26)  data,  the  A  pool  IRMMDR  time  responses  ini¬ 
tially  rise  as  AT-changing  collisions  populate  A  symmetry 
states  with  pumped  molecules.  However,  in  contrast  to  the 
9P(26)  data,  the  A  pool  time  responses  for  9/?  (12)  continue 
to  rise  as  the  pool  temperature  relaxes  downward  towards 
cell  temperature  and  the  population  increases  in  the  mea¬ 
sured  rotational  states.  Because  of  the  monotonic  rise  of 
these  A  symmetry  time  responses,  the  temperature  relaxation 
does  not  cause  a  distinctive  ETF  such  as  that  of  the  9P(26) 
data. 

This  thermal  AT-c hanging  RET  process  with  a  relaxing 
temperature  is  modeled  in  the  following  manner:  the  pro¬ 
gram  initially  calculates  the  total  population  that  AT-changing 
collisions  put  in  all  A  states  using  the  parameter  k±K ,  which 
is  the  same  for  both  pump  states.  The  IRMMDR  time  re¬ 
sponse  of  a  specific  A  symmetry  transition  is  then  calculated 
using  this  time-varying  A  pool  population  and  the  rotational 
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FIG.  9.  Measured  IRMMDR  lime  response  (dots)  of  the  J  =6-1,  K  =4,  l  =  1 
transition  at  100  mTorr.  The  solid  line  is  from  a  fit  which  includes  detector 
bandwidth  modeling,  while  the  dashed  line  is  based  on  a  model  that  neglects 
the  bandwidth  effects.  The  digitized  pump  pulse  of  the  9PI26)  line  is  shown 
at  the  bottom. 


transition’s  time-varying  absorption  coefficient.  The  tem¬ 
perature  dependence  of  a  mm/submm  absorption  coefficient 
for  CH3C1  is  given  by 

a(r)«r_5/2  exp| j,  (6) 


where  £  is  the  rotational  energy  of  the  lower  state  and  the 
limit  hv<£EkT  has  been  used  to  simplify  the  expression.45 
The  IRMMDR  time  response  of  a  specific  4  symmetry  tran¬ 
sition  is  obtained  by  multiplying  the  population  in  the  A  pool 
transported  by  the  /(-changing  rotational  process  by  the  ratio 


ait) 
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512 


exp 


£  [  1  1 

k  i  Tt  ~W) 


(7) 


where  Tr  is  the  ambient  temperature  and  T(t)  is  the  tempera¬ 
ture  of  the  population  transported  by  /(-changing  collisions. 
The  time  dependence  of  a(t)  is  obtained  by  substituting  Eq. 
(5)  into  Eq.  (7).  The  modeled  time  responses  shown  in  Fig.  6 
were  calculated  using  these  equations.  Figure  8  compares 
time  responses  of  the  unmodified  CH3F  model  (pool  tem¬ 
perature  fixed  at  ambient  temperature)  to  those  of  the  gener¬ 
alized  model  described  above.  Not  only  does  this  model  cor¬ 
rectly  demonstrate  the  /(-dependent  ETF  of  the  9P(26)  data, 
it  also  models  the  9/?  (12)  data,  which  does  not  demonstrate  a 
distinctive  ETF.  using  the  same  rate  constant  for  /(-changing 
collisions  k^K  and  the  same  relaxation  rate  constant  for  the 
temperature  kXR  . 


yields  insight  into  RET — the  population  flow  mechanisms. 
Fortunately,  most  of  the  time  responses  we  acquired  are 
dominated  by  collisional  population  flow;  however,  by  cor¬ 
rectly  modeling  the  other  two  sources,  the  quality  of  the  fits 
is  significantly  improved,  especially  when  modeling 
J -changing  collisions  within  a  few  J  of  the  pumped  state. 
One  of  the  tests  of  the  quality  of  the  model  is  whether  the 
rate  constants  are  independent  of  pressure  as  they  should  be 
for  two-body  collisions.  Correctly  modeling  the  effects  of 
Doppler  and  collisional  dephasing  as  well  as  instrumental 
bandwidth  reduced  pressure  dependencies  in  the  fits  and  thus 
allowed  us  to  apply  the  model  over  a  larger  pressure  range 
when  measuring  the  collision  rates. 

Figure  9  shows  data  significantly  affected  by  the  detec¬ 
tor’s  bandwidth.  At  these  higher  pressures,  where  the  rates 
are  fastest,  the  IRMMDR  time  responses  exhibit  peak  delay 
and  overshoot  characteristic  of  a  two-pole  lowpass  filter.  The 
detector’s  response  is  approximately 
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'  /(Oexp 
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\/4 <»l-r  , 
xstn - - -  (t-t0) 


dta. 
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where  f(t)  is  the  input  signal  measured  by  the  detector  and 
f'  =  max(20/y,  0).  We  ignore  the  constant  amplitude  factor  C 
since  amplitude  is  a  free  parameter  in  our  fits.  The  fit  ad¬ 
justed  the  oscillation  parameter  a>0  and  decay  parameter  y, 
resulting  in  values  of  0.64  and  0.59  MHz,  respectively.  Al¬ 
though  only  the  fastest  IRMMDR  time  responses  were  sig¬ 
nificantly  affected  by  the  bandwidth  (starting  at  about  20 
mTorr  for  the  states  populated  by  7-changing  collisions  and 
at  about  100  mTorr  for  those  populated  by  /(-changing  col¬ 
lisions),  by  correctly  modeling  the  detector  bandwidth,  we 
were  able  to  fit  data  from  a  larger  range  of  pressures. 

At  lower  pressures,  where  the  collisional  rates  are  slow¬ 
est,  we  observed  another  bandwidth  limiting  effect  that  is 
due  to  the  molecules'  polarization  response  times.  This  may 
be  thought  of  as  the  time  it  takes  the  induced  dipoles,  the 
pumped  nonthermal  molecules,  to  align  with  the  probe  field. 
The  response  may  be  derived  using  the  density  operator 
formalism.46  For  our  typical  probe  intensity,  the  Rabi  fre¬ 
quency  is  8000/271  Hz,  which  is  slow  compared  to  the  time 
scales  of  our  measurements  and  is  much  slower  than  the 
dephasing  rates;  therefore,  the  calculation  is  performed  in  the 
limit  that  the  Rabi  frequency  is  much  longer  than  the  relax¬ 
ation  rates.  We  want  to  find  the  time-dependent  absorption 
a(t)  of  the  probe  field 


dP 

dz 


~ a(t)P , 
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B.  Other  model  modifications 

Several  sources  contribute  to  the  probe  time  responses 
that  we  measure:  (1)  population  changes  in  either  the  upper 
or  lower  rotational  level;  (2)  dephasing  due  to  Doppler 
broadening  and  pressure  broadening  from  both  elastic  and 
inelastic  collisions;  and  (3)  instrumental  bandwidth.  Thus  far, 
we  have  been  concerned  with  the  source  that  most  directly 


where  z  is  the  direction  of  propagation  of  the  radiation  along 
the  axis  of  the  absorption  cell  and  P  is  the  mm/submm  probe 
power.  The  electric  field  in  the  molecular  frame  is 

=  +  c.c.  =  T%(t)e~'iu,~iv)‘  +  c.c. 

(10) 


and  the  polarization  is 
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FIG.  10.  Measured  IRMMDR  time  response  (dots)  of  the  7=5-6,  K= 4, 
/  =  I  transition  at  20  mTorr.  The  solid  line  is  based  on  the  complete  model 
which  includes  the  modeling  of  polarization  response  time,  while  the  dashed 
line  is  based  on  a  model  which  neglects  polarization  delay.  The  digitized 
pump  pulse  of  (he  9/>(26)  line  is  shown  at  the  bottom. 


P(v,t)  =  ^(v,t)e-,(,u-tv)'  +  c.c.,  (11) 

where  %(t)  and  ^/‘(v,t)  are  the  slowly  varying  amplitudes 
in  the  molecule  frame.  The  power  per  unit  volume  absorbed 
from  the  probe  field  as  measured  by  the  detector  is 

dP  dl  f 

dV=Z^J  dv(f(v,t)-E(v,l)), 

=  -jf  dv  Rt[-i(o%*(t)  sAv,t)],  (12) 

where  V  is  the  cell  volume  and  /  is  the  probe  intensity.  The 
velocity  integration  is  performed  over  the  probed  velocity 
subclass.  The  polarization  may  be  derived  using  the  equation 
of  motion  for  the  density  matrix  p, 

p(i/,f)=--  [H,p(v,t)],  (13) 

where  and  the  density  operator  is  a  func¬ 

tion  of  molecular  velocity  v  and  time.  Labeling  the  lower 
and  upper  probe  states  1  and  2,  respectively,  we  expand  Eq. 
( 13)  and  modify  the  result  to  include  decay  from  dephasing 
collisions 

i 

P2l(V,t)=  -(ypB.  +  t'a»2l)P2l+£  J*2i-E(l7,f)[Pu(t;,f) 

~Pn(v,t)],  (14) 

pl2(v,t)  =  P2i(V,t) 

where  is  the  mm/submm  transition  frequency  between 
the  rotational  levels  and  >B  is  the  half-width  at  half¬ 
maximum  (HWHM)  pressure  broadening  rate  in  radians  per 
seconds.  The  solution  to  Eq.  (14)  in  the  rotating  wave  ap¬ 
proximation  is 


7*21  -e£n 

p2i(v,t)  =  i  —  I  ^  dt'  exp[  —  (  ypB  +  'u»2i)(/-/')] 

Xexp[  —  t(  o>  — Art)  )f'][p  |  ,(o,r')  -p22(u,f')], 

,  .  *  .  .  (15) 

p12(t;,r)  =  p21(i;,f). 

The  polarization  is  given  by 

P(u,/)  =  Tr[/ip(t;,/)]  =  2  Re[/a,2P2i(u,/)]-  (16) 

Compared  to  Eq.  (11),  the  slowly  varying  amplitude  of  the 
polarization  can  be  written  as 

yAv,t)  =  2pLnp2l(v,t)e‘{u~kuU.  (17) 

Rearranging  Eq.  (9)  and  substituting  Eqs.  (12),  (15),  and  (17) 
yields  the  general  expression  for  the  time  dependence  of  the 
absorption  coefficient 

4 itA:|/*2i  -e|2  /  f  f> 

a(t)  = - h -  Re(J  dvj. 

x  exp{ -[yPB.+i(  A -A: t))](/-r')}[p, ,(!;,/') 

(18) 

where  A=o)12-oj  and  Ar=ti>c~’.  For  A=0  and  assuming 
a  thermal  distribution  of  molecular  velocities,  i.e., 
pu(v,t)  =  pu(t)(u\fn)~'  exp(-i;2/u2),  this  reduces  to 

47rit|/tt2i -e|2  ( [i 

a(,)  = - h  Re  J  dt'e-v* * 


/ 1  \2  ,1 

Xexp 

-  2  *«  ('~f') 

[P^(,,)-P22(,,)] 

(19) 

where  u  =  sj2kT/M  is  the  thermal  speed  of  the  molecules 
along  the  axis  of  the  absorption  cell.  The  effect  that  the  po¬ 
larization  response  has  on  a  pumped  K/l  stack  IRMMDR 
time  response  is  shown  in  Fig.  10.  This  model  is  not  com¬ 
plete  since  it  assumes  a  thermal  velocity  distribution  for  mol¬ 
ecules;  nonetheless,  including  this  polarization  effect  in  the 
model  significantly  improved  the  quality  of  the  fits  without 
adding  any  adjustable  parameters.  Data  significantly  affected 
by  nonthermal  velocity  distributions  were  excluded  from  the 
fits,  specifically,  low  pressure  J-changing  RET  data  near  the 
pumped  state.  Allowing  the  model  to  adjust  the  pressure 
broadening  parameter  in  Eq.  (19)  resulted  in  the  value 
ypB.= 15  MHz/Torr,  which  agrees  well  with  pressure  broad¬ 
ening  coefficients  based  on  linewidth  measurements;  there¬ 
fore,  ypB.  was  fixed  ,0  this  value. 

In  order  to  ensure  that  no  pertinent  early  time  informa¬ 
tion  was  inadvertently  lost  as  a  result  of  undersampling,  we 
typically  used  digitization  rates  of  0.133  and  0.333  p s  when 
acquiring  IRMMDR  time  responses.  When  analyzing  the 
data,  we  employed  a  scheme  to  eliminate  oversampied  data 
points  from  the  fit  to  reduce  computer  run  times.  This  was 
accomplished  by  programming  the  computer  model  to  apply 
a  boxcar  averager  with  a  window  size  depending  on  the  fre¬ 
quency  content  of  the  data  trace  at  a  given  time.  Problems 
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with  changing  the  phase  of  the  signals  were  avoided  by  ap¬ 
plying  the  boxcar  averager  to  the  measured  time  responses 
after  simulated  time  responses  had  already  been  subtracted 
from  them  in  preparation  for  the  least  squares  fit  routine.  The 
windows  were  required  to  be  narrow  enough  that  the  filter 
affected  fitted  parameters  by  no  more  than  one  part  in  10'. 
Even  with  this  stringent  constraint,  many  fit  times  were  re¬ 
duced  by  more  than  90%. 

C.  Application  of  the  model 

Excluding  data  acquired  with  various  pump  frequency 
and  probe  frequency  offsets  (used  for  observing  velocity  sub¬ 
class  effects),  there  were  1027  IRMMDR  time  responses 
available  for  the  CH3  3SC1  analysis.  Where  redundant  time 
responses  existed  (same  pump  line,  probe  line,  and  pressure), 
the  one  with  the  best  signal  to  noise  ratio  was  chosen,  leav¬ 
ing  579  time  responses  for  64  probe  transitions  and  two 
pump  coincidences.  Most  probe  transitions  were  measured 
for  both  pump  states,  and  most  time  responses  were  taken  at 
six  to  eight  different  pressures.  After  eliminating  the  IRM¬ 
MDR  time  responses  with  inadequate  signal  to  noise  ratios 
and  many  of  the  time  responses  for  vibrational  processes, 
there  were  366  time  responses  that  were  weighted  in  the  fits 
according  to  their  signal  to  noise  levels. 

To  determine  the  collision  rates  for  CH3C1,  IRMMDR 
time  responses  were  grouped  by  pool  or  included  as  indi¬ 
vidual  7  states  for  the  pumped  K/l  stack,  then  the  following 
sequence  applied:  (1)  a  large  set  of  time  responses  over  the 
full  range  of  pressures  (up  to  297  time  responses  simulta¬ 
neously)  was  fit  by  adjusting  the  appropriate  parameters. 
Eventually,  all  data  were  grouped  into  two  sets  (one  for  all 
time  responses  of  states  in  the  same  K/l  stacks  as  the 
pumped  states  and  one  for  all  other  data).  Since  the  program 
runs  much  faster  when  the  individual  7  states  are  not  in¬ 
cluded,  the  model  used  individual  7  states  only  for 
J  -changing  data  and  the  faster  model  using  only  pools  of 
states  was  used  for  all  other  time  responses.  (2)  A  data  set 
was  subdivided  into  groups  with  the  same  pressure  and  then 
were  fit  separately  to  check  for  any  pressure  dependence  in 
the  parameters.  (3)  Since  pressure  variations  in  the  fitted  rate 
constants  indicate  a  systematic  experimental  or  modeling 
problem,  data  for  any  fits  demonstrating  deviations  with 
pressure  at  either  low  or  high  pressures  are  eliminated.  How¬ 
ever.  the  model  was  found  to  be  valid  over  a  large  pressure 
range  and  only  16  time  responses  for  some  transitions  in  the 
pumped  K/l  stack  at  low  pressures  were  eliminated  due  to 
significant  collisionless  effects.  (4)  Finally,  the  fit  was  run 
again  using  the  entire  data  set  to  arrive  at  the  final  param¬ 
eters. 

These  steps  were  performed  for  each  group  of  data  and 
then  the  whole  process  was  repeated  until  the  fit  error  was 
minimized.  The  final  fit  parameters  represent  350  IRMMDR 
time  responses  grouped  in  two  sets,  each  of  which  were  fit 
simultaneously— 53  time  responses  representing  7-changing 
RET  and  297  time  responses  representing  K -changing  RET 
and  the  vibrational  processes.  For  K -changing  RET  and  vi¬ 
brational  processes,  the  model  was  found  to  be  valid  over  the 
full  pressure  range  used  to  acquire  data.  1-100  mTorr  for 
K -changing  RET  and  20-1000  mTorr  vibrational  processes. 


time  (|is) 


FIG.  11.  Measured  IRMMDR  time  response  (dotsl  of  the  u6=l,  7=7-8, 
K=  4,  f  =  l  transition  at  several  pressures.  The  solid  lines  are  based  on  the 
complete  model  using  the  IOS-P  law  for  7 -changing  collisions.  The  digi¬ 
tized  pump  pulse  of  the  9PI26)  line  is  shown  at  the  bottom. 


For  7 -changing  RET,  the  dipole-dipole  process  and  the 
scaling  laws  were  applied  to  the  IRMMDR  time  responses  of 
the  pumped  K/l  stacks  to  determine  the  best  method  of  mod¬ 
eling  J -changing  collisions.  Four  permutations  of  the  model 
were  run  using  three  groupings  of  the  data  [only  9P(26) 
data,  only  9R(12)  data,  and  both  simultaneously]  to  deter¬ 
mine  the  best  model— the  SPG  law  with  and  without  a  sepa¬ 
rate  dipole-dipole  process  and  similarly  the  IOS-P  law  with 
and  without  a  separate  dipole-dipole  process.  When  the 
dipole-dipole  process  was  included  in  the  model,  the  scaling 
law  was  used  for  only  higher-order  collisions  A7^2,  other¬ 
wise,  the  scaling  law  included  A735 1 .  The  parameters  for  the 
bandwidth  modeling  and  the  polarization  response  were  also 
adjusted  in  these  A7  fits  and  the  resulting  parameters  were 
then  used  for  all  subsequent  fits.  We  came  to  several  conclu¬ 
sions.  First,  both  the  SPG  and  IOS-P  scaling  laws  accurately 
model  the  data  with  the  latter  resulting  in  a  slightly  better 
figure  of  merit.  Also,  we  found  that  fitting  separately  for  the 
dipole-dipole  process  did  not  significantly  improve  the  qual- 


TABLE  11.  RET  fit  parameters. 


Parameter 


Fit  result' 


Cross  section  (A2) 


‘as 

11. 110.8  us  ’TorT1 

ToMPtU) 

13816  K 

C..wui 

■  374117  K 

^T» 

7.510.2  jisMTorr'' 

4812  jis''TorT'’ 

IOS-P  c 

21.711.1  jjs-'Torr'1 

y 

1.2310.02 

SPG  c 

0.012010.0006  ^u  'Torr' 

r 

1.1210.01 

6915 


4711 
299  ±9 
13417 


0.07410.005 


'Errors  are  2.5  times  <r  from  the  nonlinear  least  squares  fit  routine. 
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iiy  of  the  fits  and  is  therefore  unnecessary  to  include  in  the 
model.  That  is.  either  the  SPG  or  10S-P  scaling  law  accu¬ 
rately  accounts  for  the  A7  =  ±l  collisions  in  addition  to  the 
higher-order  |A7|^2  collisions.  Examples  of  the  fits  using 
the  IOS-P  model  are  shown  in  Figs.  5  and  6  for  different  7 
states  and  in  Fig.  11  at  different  pressures.  The  total  rate 
constants  for  dipole-dipole  collisions  using  the  branching 
ratios  of  Eq.  (1)  is  /tdd =48  /zs"1  Torr-1,  about  six  times  the 
gas  kinetic  collision  rate  constant. 

In  the  previously  reported  analysis  of  the  CH3F  data.6 
the  dipole-dipole  process  was  included  separately.  However, 
with  the  improved  bandwidth  modeling  and  the  new  polar¬ 
ization  response  modeling  in  the  revised  model  discussed 
here,  a  reanalysis  of  the  data  used  in  the  previous  paper 
indicates  that  a  separate  dipole-dipole  process  is  unneces¬ 
sary  for  CH3F  as  well.  These  findings  will  be  reported  in 
detail  in  a  future  paper. 

The  resulting  RET  parameters  for  both  7-  and 
AT-changtng  collisions  are  listed  in  Table  II.  Rates  are  con¬ 
verted  to  cross  sections  using  the  relation 

k  =  n(v,el)tx, 

where  n  is  the  number  of  molecules  per  volume  and  (ure|)  is 
the  average  relative  velocity  of  the  collision  partners.  The 
average  relative  velocity  for  the  Maxwell-Boltzmann  distri¬ 
bution  is 

<«/„,)=  Vl6*77  7 rm. 

All  parameters  come  from  fits  simultaneously  including  data 
from  both  pump  coincidences.  The  scaling  law  parameters 
for  both  the  SPG  and  IOS-P  models  come  from  fits  where  the 
selected  scaling  law  modeled  all  7-changing  processes  in¬ 
cluding  those  with  dipole-dipole  selection  rules.  The  kM  pa¬ 
rameter  comes  from  a  fit  where  the  IOS-P  law  modeled  all 
other  7-changing  collisions,  specifically  |A7|st2,  while  Eq. 
(1)  was  used  for  A7  =  ±l. 

D.  Other  possible  sources  of  A K  signal  considered 

Since  a  central  result  of  this  study  is  the  generalization 
of  the  thermal  K -changing  rotational  process  to  accurately 
model  the  data  for  both  pump  coincidences,  it  is  important  to 
consider  other  possible  sources  of  the  ETF.  One  such  source 
is  velocity  subclass  effects.  The  center  frequencies  of  the 
pump  transitions  do  not  exactly  match  the  center  frequencies 
of  the  pumped  CH3C1  rovibrational  transitions — the  laser  fre¬ 
quency  mismatches  are  20  MHz  for  the  9P(26)  coincidence 
with  CH335C1  and  -30  MHz  for  9R(12).47  Since  the  CH3C1 
infrared  transitions  are  inhomogeneously  broadened  in  our 
pressure  range,  the  molecules  that  are  pumped  are  those  with 
the  velocity  that  provides  the  needed  Doppler  shift  to  correct 
for  this  frequency  mismatch.  The  pumped  molecules  retain 
this  nonthermal  velocity  distribution  until  thermaiized  by  gas 
kinetic  collisions.  Since  the  dipole-dipole  and  cumulative 
rate  for  K -changing  collisions  are  comparable  to  or  faster 
than  the  gas  kinetic  collision  rate,  many  molecules  undergo  a 
change  in  rotational  state  before  their  velocity  distribution  is 
thermaiized.  So  the  frequency  spike  that  is  observable  in  the 
probe  frequency  profiles  of  the  three  level  double  resonance 
measurements  may  also  be  observed  for  other  probe  transi¬ 


tions,  specifically,  transitions  in  the  pumped  K/l  stack  and  to 
a  lesser  degree  all  transitions  of  the  same  symmetry  as  the 
pumped  state.  In  their  infrared/infrared  double  resonance 
study  of  CH3F,  Song  and  Schwendeman  observed  these 
spikes  in  the  frequency  profiles  of  IRIRDR  measurements  of 
the  pumped  K  stack.15  Everitt  also  observed  this  effect  for 
both  7-  and  K -changing  collisions  in  time-resolved 
IRMMDR  measurements.37  To  investigate  the  possibility  that 
the  ETF  of  the  CH3C1  data  is  due  to  this  velocity  subclass 
effect  IRMMDR  time  responses  were  taken  over  a  grid  of 
pump  and  probe  frequencies.  For  our  mm/submm  probe  fre¬ 
quencies,  it  is  only  at  pressures  less  than  about  10  mTorr  that 
pressure  broadening  becomes  smaller  than  Doppler  broaden¬ 
ing.  Although  there  is  marginal  evidence  for  an  early-time 
probe  frequency  asymmetry  in  the  A  transitions  at  lower 
pressures,  the  size  of  the  effect  is  clearly  smaller  than  the 
ETF  of  Fig.  6.  Furthermore,  the  fact  that  the  magnitude  of 
the  ETF  increases  at  higher  pressures  is  not  consistent  with 
attributing  the  ETF  to  velocity  effects.  Moreover,  the  absence 
of  the  ETF  in  the  9/?  (12)  data  is  not  consistent  with  attrib¬ 
uting  the  ETF  to  velocity  subclass  effects. 

We  also  attempted  to  account  for  the  ETF  of  the  9P(26) 
IRMMDR  time  responses  by  breaking  K -changing  RET  into 
separate  processes  K'  =  3— -0,  K'=  3— >6,  and  K'  =  0— *6.  As 
with  the  CH3F  model,  the  individual  A K'  processes  were 
defined  to  be  thermal  in  7  with  the  temperature  fixed  at  the 
cell  temperature.  The  reasoning  behind  the  model  was  that 
the  ETF  for  the  K'  =  0  and  K'  =  3  data  might  be  caused  by 
population  initially  entering  these  states  quickly  then  moving 
into  K'=s 6  states  more  slowly.  Although  such  a  model  can 
qualitatively  account  for  some  of  the  features  observed  in  the 
data,  it  could  not  reproduce  the  IRMMDR  time  responses 
due  to  the  following  competing  effects:  the  individual  A K' 
rates  would  have  to  be  close  enough  to  each  other  to  account 
for  the  similar  rise  times  of  K'=  0,  3,  and  6  transitions,  but  at 
the  same  time  different  enough  to  allow  for  the  fall  of  the 
ETF  in  the  /C'=0  and  3  data  that  is  absent  in  the  K' =6  data. 
Given  the  equilibrium  percentage  of  population  in  K'~>6 
relative  to  all  u6=  1,  A  symmetry  states  (31%),  no  combina¬ 
tion  of  rates  could  account  for  the  ETF  in  the  K'  =0  and  3 
states.  Further  subdividing  the  A K'  process  into  more  state- 
specific  processes  also  failed  to  overcome  this  difficulty. 

VII.  DISCUSSION 

A.  Comparison  of  CH3CI  and  CH3F 

Within  the  expected  rate  constant  variations,  the  rota¬ 
tional  processes  in  CH3CI  are  for  the  large  part  the  same  as 
those  of  CH3F.  The  dipole-dipole  process  has  the  fastest 
state-to-state  rates  in  both  molecules;  however,  other  RET 
processes  in  both  molecules  play  more  important  roles  in  the 
thermalization  of  the  rotational  states  in  the  pumped  vibra¬ 
tional  manifold  since  they  have  larger  ranges  in  A7  and  per¬ 
mit  K -changing  collisions.  Higher-order  7-changing  colli¬ 
sions  in  both  molecules  populate  states  in  the  pumped  K 
stack  (or  K/l  stack  for  p6=l  in  CH3C1)  that  are  more  than 
several  7  away  from  the  pumped  state  more  quickly  than 
multiple  dipole-dipole  collisions  because  they  have  a  larger 
range  in  7.  In  addition,  the  thermal  K -changing  RET  process 
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in  both  molecules  exhibits  state-to-state  rates  very  different 
in  character  than  those  of  /-changing  collisions.  Although 
some  of  the  double  resonance  signals  for  CH3C1  and  CH3F 
appear  to  be  quite  different,  this  disparity  can  be  explained 
by  a  difference  in  the  temperatures  at  which  X -changing  col¬ 
lisions  populate  states  of  the  pumped  symmetry  in  the  two 
molecules. 

Due  to  the  similarity  of  the  molecules,  the  generalized 
model  developed  for  CH,C1  should  also  apply  to  CH3F.  In 
fact,  a  preliminary  reanalysis  of  the  C H3F  data  indicates  that 
the  generalized  model  is  consistent  with  CH3F.  However,  the 
previous  work5  found  that  a  model  using  a  fixed  300  K  tem¬ 
perature  for  the  thermal  Boltzmann  distribution  of  final  states 
was  sufficient  for  modeling  X -changing  collisions  in  CH3F, 
whereas  fixed  temperature  models  fail  for  CH3C1.  Several 
reasons  for  this  follow:  First,  since  the  B  rotational  constant 
for  CH3C1  is  nearly  half  the  value  of  B  for  CH3F,  while  their 
A  rotational  constants  are  nearly  equal,  there  are  many  more 
thermally  populated  7  states  for  a  given  X  in  CH3C1  than  in 
CH3F.  Consequently,  CH3C1  provides  a  more  severe  test  of 
the  thermal  nature  of  X -changing  RET  since  large  changes  in 
the  angular  momentum  of  the  pumped  population  are  re¬ 
quired  to  warm  or  cool  the  population  distribution  of  A  pool 
states  to  room  temperature.  Second,  the  available 
C02-CH3CI  pump  coincidences  populate  rovibrational  states 
that  are  more  widely  spaced  in  rotational  energy  than  those 
of  CH3F.  For  CH3C1,  the  pumped  rovibrational  states  have 
rotational  energies  of  85  cm"1  for  9P{26)  and  285  cm-1  for 
9/?(12),  while  for  CH,F,  the  pumped  rovibrational  states 
have  rotational  energies  of  64  cm"1  for  9P(32)  and  148 
cm"1  for  9P(20).  This  difference  of  pumped  state  rotational 
energies  in  CH3C1  and  the  fact  that  one  is  larger  than  room 
temperature  thermal  energy  while  the  other  is  smaller  en¬ 
hances  the  differences  the  molecules  exhibit  in  X -changing 
collisions.  Last,  a  secondary  pump  coincidence  for  the 
9<P(20)  pump  line  noticeably  affects  early  times  in  IRM- 
MDR  time  responses  needed  for  measuring  X -changing  col¬ 
lisions  in  AT  =  1  of  12CH3F. 

B.  Comparison  with  NH3 

In  addition  to  the  CH3F  work5  and  the  present  study,  the 
work  of  Abel  et  al.  for  NH3  provides  detailed  measurements 
of  state-to-state  rates  for  X-changing  collisions  for  a  sym¬ 
metric  top  molecule.29  Using  a  time-resolved  IRIRDR  tech¬ 
nique,  they  probed  the  time  responses  of  several  rotationally 
resolved  2v2«—  transitions  in  NH3  as  the  states  responded 
to  excitation  by  the  10P(32)  line  of  a  C02  laser.  This  pump 
coincidence  populates  the  a  inversion  level  of  u2=l,  7=5, 
X  = 3.  They  probed  several  rotational  states  of  opposite 
parity — six  7  levels  in  X  =  3  and  eight  /  levels  in  X  = 0.  This 
allowed  them  to  calculate  the  state-to-state  RET  rate  con¬ 
stants  for  several  AX=0  and  -3  collisional  processes  using 
a  master-equation  model  for  RET.  In  their  model,  the  AX' =0 
processes  were  calculated  using  the  EGL  fitting  law,  but  the 
A/C  =  -3  processes  were  fitted  individually  since  neither  the 
EGL  nor  SPG  law  could  reproduce  satisfactorily  the  IRIRDR 
time  responses.  They  grouped  many  rotational  states  with 
larger  7  or  X  into  four  pools  of  states  to  reduce  the  total 
number  of  states  in  their  model  since,  in  the  absence  of  a 


(a)  CH3CI 


(b)  NH3 


FIG.  12.  Slate-to-stale  collision  cross  sections  for  rotational  collisions  in  (a) 
CHjCI  using  the  9PI26)  pump  line  and  (b)  parity-changing  rotational  colli¬ 
sions  measured  by  Abel  el  al.  (Ref.  29)  for  NH,.  The  cross  sections  are 
normalized  to  each  molecule’s  gas  kinetic  cross  section.  The  initial  state  for 
CH,Cl  is  u*=l,  J- 6,  K= 4,  f=l  and  the  initial  state  for  NHj  is  the  a 
inversion  level  of  y,=  1,  J  =  5.  K  =  3. 


scaling  relation  for  X-changing  collision  rates,  each  addi¬ 
tional  state  adds  a  free  parameter  to  the  master  equation. 

It  is  interesting  to  compare  the  magnitudes  of  the  state- 
to-state  collision  rate  constants  for  CH3C1  and  NH3,  bearing 
in  mind  that  while  both  molecules  are  symmetric  tops,  NH3 
has  the  additional  degree  of  freedom  from  inversion.  Abel 
et  al.  found  that  both  parity-changing  (inversion  state¬ 
changing)  and  parity-conserving  (inversion  state-conserving) 
collisions  occurred  and  modeled  them  by  assuming  their 
rates  are  different  by  a  ratio  0=k(a>->a)/k(a+-*s ) 
=  k(s'—s)/k(a*-*s),  which  was  found  to  be  greater  than  0.7. 
Figure  12  compares  the  state-to-state  rate  constants  for 
parity-changing  collisions  in  NH3  measured  by  Abel  et  al 
with  the  state-to-state  rate  constants  for  collisions  in  CH3G 
calculated  using  our  five  parameter  model  for  RET.  The  rates 
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have  been  normalized  to  the  gas  kinetic  collision  cross  sec¬ 
tions  for  each  molecule. 

Remembering  that  7,  =7  pump=6  for  CH3C1  and 
7,  =7pump=5  for  NHj,  Fig.  12  shows  that  for  both  molecules, 
the  largest  state-to-state  rates  are  those  with  dipole-dipole 
selection  rules.  In  addition,  for  NH3,  the  fastest  A/(=  -3 
state-to-state  rate  is  A7 =0  and  the  rates  decrease  with  larger 
|A/|,  indicating  a  propensity  for  RET  processes  to  conserve 
angular  momentum  in  the  pumped  molecule  during 
K -changing  collisions.  In  contrast,  the  state-to-state  rates  for 
CH3C1  AT -changing  collisions  are  maximum  at  7  =  10  rather 
than  7  =6  because  for  this  molecule  /(-changing  RET  is  ther¬ 
mal;  7  =  10  corresponds  to  the  maximum  in  the  Boltzmann 
distribution  for  the  138  K  initial  temperature  of  the  9/*(26) 
pump  coincidence  shown  in  Table  II.  Final  states  at  smaller 
Jf  receive  less  population  because  their  (27+ 1)  M - 
degeneracies  are  small,  while  states  at  larger  7 f  receive  less 
population  due  to  the  exponential  rolloff  of  the  Boltzmann 
distribution  with  energy. 

Figure  12  shows  that  the  7-changing  collisions  in  CH3C1 
are  faster  than  the  fastest  /(-changing  collisions  by  a  factor 
of  about  10  (for  large  A7)  to  100  (for  small  A7),  whereas  in 
NH3,  A/(=0  collisions  are  only  two  to  three  times  faster 
than  A/(=— 3  collisions.  These  ratios  can  be  understood 
from  the  physical  structures  of  the  two  molecules  using  the 
following  classical  picture:  collisions  that  change  K  for  both 
CH3C1  and  NH3  require  a  collision  of  hydrogen  atoms  since 
they  provide  the  off-axis  structure  for  these  symmetric  tops. 
The  center  of  mass  for  NH3  is  near  the  nitrogen  nucleus,  and 
the  rotational  constants  for  NH3  are  large  and  similar  in  value 
"»300  GHz  and  C=>190  GHz);  thus,  the  moment  of  inertia 
for  rotation  about  the  symmetry  axis  is  not  significantly  dif¬ 
ferent  than  the  moment  of  inertia  about  a  perpendicular  axis 
for  tumbling  motion.  While  7-changing  collisions  can  occur 
for  long-range,  dipole-dipole  interactions  of  NH3  molecules, 
/(-changing  collisions  require  short-range,  higher-order  inter¬ 
actions  closer  to  the  repulsive  portion  of  the  intermolecular 
potential.  During  such  a  hard  collision  between  the  partner 
and  a  hydrogen  atom  of  NH3,  it  is  about  equally  likely  that 
the  molecule’s  angular  momentum  about  the  symmetry  axis 
will  change  as  for  an  axis  of  tumbling  motion  because  the 
line  defined  by  the  point  of  contact  and  the  molecule’s  center 
of  mass  approximately  bisects  the  angle  between  these  two 
rotation  axes.  Therefore,  while  A/(=0  state-to-state  rates  can 
be  expected  to  be  larger  since  the  dipole -dipole  interaction 
between  the  colliding  molecules  do  not  require  small  impact 
parameters,  during  short-range  collision  of  NH3  molecules 
where  higher-order  multipoles  can  be  significant,  7-  and  K- 
changing  collisions  are  roughly  equally  likely  for  NH3  as 
observed  in  Fig.  12. 

The  structure  of  CH3C1  is  different  than  that  of  NH3 .  The 
center  of  mass  for  CH3G  is  near  the  heavy  chlorine  nucleus 
and  the  moments  of  inertia  are  different  by  more  than  a  fac¬ 
tor  of  10  13  GHz  and  A  *>150  GHz).  So  a  hard  collision 

between  a  partner  and  a  hydrogen  atom  in  a  CH3G  molecule, 
the  type  of  collision  required  to  change  K,  is  more  likely  to 
change  the  angular  momentum  about  an  axis  for  tumbling 
motion  than  about  the  symmetry  axis  since  the  molecule's 
center  of  mass  is  at  the  opposite  end  of  the  molecule.  In  other 


pressure  (mT) 


FIG.  13.  Measured  pressure  broadening  linewidths  vs  pressure  for  the 
c/s=l,  7=6-7,  K= 4.  /  =  1  transition.  The  inset  shows  one  of  the  frequency 
profiles  taken  at  20  mTorr.  The  Voigl  line  shape  was  used  with  the  Gaussian 
component  fixed  at  0.163  MHz.  Both  the  intercept  and  slope  were  free 
variables  in  the  linear  fit,  resulting  in  a  pressure  broadening  coefficient  of 
14.3±1.0  MHz/Torr. 


words,  a  CH3C1  collision  that  changes  K  is  likely  to  have 
caused  an  even  larger  change  in  7.  In  fact,  such  large 
changes  in  7  occur  that  the  /(-changing  process  in  CH3G  can 
be  treated  as  a  thermal  process. 


C.  A  comparison  of  collision  rates  with  pressure 
broadening 

An  important  cross  check  of  the  accuracy  of  the  mea¬ 
sured  state-to-state  collision  rates  is  a  comparison  of  the 
pressure  broadening  derived  from  the  individual  rate  con¬ 
stants  with  the  measured  pressure  broadening  coefficients. 
For  CH3G,  the  cumulative  depopulation  rate  of  the  9P(26) 
pumped  state  through  AK=0  channels  is  85  /us-1  Torr-1 
based  on  an  IOS-P  law  summation.  This  includes  the  A7=0, 
A/(=0  elastic  collision  rate  that  is  calculated  using  the 
branching  ratios  for  dipole-dipole  collisions  found  in  Eq. 
(1).  It  is  interesting  to  note  the  total  rate  constant  for  dipole- 
dipole  collisions  is  48  Ton--1,  only  57%  of  all  AK=0 

collisions.  As  shown  in  Table  II,  the  cumulative  rate  constant 
for  /(-changing  collisions  is  11.1  1  Ton-1.  The  equation 

used  to  calculate  pressure  broadening  in  MHz/Torr  is 

y/p=(  1/2  *■)(*„+*, )/2,  (20) 

where  p  is  the  pressure  in  Torr  and  the  cumulative  depopu¬ 
lation  rate  constants  for  the  upper  and  lower  transition  levels 
ku  and  k,  are  expressed  in  ps~x  Torr-1.  Including  the  small 
contributions  from  the  vibrational  processes,  the  fit  param¬ 
eters  result  in  a  pressure  broadening  coefficient  of  12.9  ±0.5 
MHz/Torr.  We  have  performed  a  series  of  pressure  broaden¬ 
ing  measurements  on  a  transition  involving  the  pumped  state 
o6=l,  7=6— >7,  K =4,  1=1.  The  data  are  plotted  in  Fig.  13 
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and  the  resulting  pressure  broadening  parameter  of  14.3  ±  1.0 
MHz/Torr  agrees  well  with  the  parameter  determined  from 
the  rate  constants. 

VIII.  CONCLUSIONS 

To  summarize,  we  report  that  all  state-to-state  rotational 
energy  transfer  processes  in  vb=  1  of  CH3CI  can  be  described 
in  terms  of  a  small  number  of  physically  meaningful  param¬ 
eters.  The  CH3C1  state-to-state  rotational  collision  rates  from 
either  of  the  two  pumped  states  to  any  other  A  symmetry 
rotational  state  in  uA=l  may  be  calculated  using  the  RET 
model  described  in  this  paper  and  measured  parameters  of 
Table  11.  For  each  pump  state,  only  five  parameters  are 
needed  to  model  all  of  these  state-to-state  rates — two  param¬ 
eters  for  either  the  SPG  or  1GS-P  scaling  law  for  7-changing 
RET  and  three  parameters  for  the  temperature-relaxing,  ther¬ 
mal  AT -changing  process.  Furthermore,  four  of  these  param¬ 
eters  are  the  same  for  all  initial  states:  only  the  initial  tem¬ 
perature  of  the  thermal  A-changing  process  varies  with 
initial  state. 

It  should  be  emphasized  that  these  parameters  account 
correctly  for  a  large  number  of  experimental  double  reso¬ 
nance  time  responses  (which  often  have  very  different 
shapes)  taken  over  a  wide  range  of  parameter  space.  Specifi¬ 
cally,  350  IRMMDR  time  responses  were  used  with  pres¬ 
sures  ranging  from  1  to  100  mTorr  (20  to  1000  mTorr  for 
states  populated  by  vibrational  processes)  for  43  probe  tran¬ 
sitions  covering  a  wide  range  of  J  and  A,  pumped  and  un¬ 
pumped  symmetry  species,  and  two  different  pump  states.  As 
for  our  earlier  CH3F  (Refs.  5  and  6)  study,  the  statistical 
power  gap  fitting  law  and  the  infinite  order  sudden  approxi¬ 
mation  using  a  power  law  expression  for  the  basis  rates  ac¬ 
curately  model  J -changing  t. A/C  =0)  collisions.  Furthermore, 
the  intermolecular  potentials  for  these  two  methyl  halides 
result  in  A' -changing  collision  rates  that  can  be  modeled  very 
simply.  The  state-to-state  rates  to  all  eligible  final  states  (i.e., 
states  of  the  same  symmetry  and  vibrational  energy  as  the 
pumped  state)  are  such  that  the  pumped  population  trans¬ 
ferred  to  these  states  is  distributed  according  to  a  Boltzmann 
distribution  even  at  the  earliest  times.  The  most  significant 
difference  between  rotational  energy  transfer  in  these  two 
molecules  is  that  whereas  a  temperature  fixed  at  the  ambient 
temperature  of  the  cell  was  found  to  be  sufficient  for  model¬ 
ing  A-changing  collisions  in  CH3F,  the  CH3C1  data  required 
a  time-varying  temperature.  However,  both  parameters  for 
this  temperature  are  closely  related  to  other  known  physical 
characteristics  of  the  molecule — the  initial  temperature  is  re¬ 
lated  to  the  rotational  energy  of  the  pumped  state  and  the  rate 
of  the  temperature’s  exponential  relaxation  back  to  the  am¬ 
bient  temperature  of  the  cell  is  slightly  less  than  the  cumu¬ 
lative  rate  of  A-changing  collisions. 
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